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LONG-TERM  BIOEFFECTS  OF  435-MHz  RADIOFREQUENCY  RADIATION 
ON  SELECTED  BLOOO-BORNE  ENDPOINTS  IN  CANNULATED  RATS 
Volume  4.  Plasma  Catecholamines 

I.  INTRODUCTION 

During  the  past  50  years#  the  United  States  has  witnessed  a  period  of 
explosive  growth  In  the  radar  and  communications  fields.  This  growth  has 
Increased  the  demand  for  available  bandwidth  and  thus  has  pushed  radar  and 
communications  frequencies  Into  higher  and  higher  ranges.  Higher  frequency 
ranges  have  permitted  faster  data  transmission  rates  and  reduced  Intersystem 
electromagnetic  Interference.  However#  these  advances  have  come  at  the  expense 
of  altering  the  planet's  radiofrequency  radiation  (RFR)  environment.  Until  the 
advent  of  advanced  radar  and  communications#  cosmic  rays  and  background 
radiation  were  the  primary  sources  of  the  Earth's  electromagnetic  environment. 
Radar  and  communications  transmissions  have  since  Increased  the  electromagnetic 
background  or  ambient  radiation  at  the  planet's  surface  by  several  orders  of 
magnitude.  At  this  time#  the  biological  effects  of  exposure  to  this  omnipresent 
electromagnetic  environment  are  not  well  understood#  despite  studies  conducted 
over  the  past  several  decades. 

This  report  presents  the  results  of  plasma  catecholamine  (norepinephrine# 
epinephrine#  and  dopamine)  assays  of  blood  samples  drawn  from  a  large  population 
of  male  Sprague-Dawley  rats  exposed  to  a  1.0  mW/cm^#  435-MHz  pulsed-wave  (1.0  ps 
pulse  width#  1-kHz  pulse  rate)  RFR  environment  for  a  6-month  duration.  The 
exposure  group  consisted  of  100  cannulated  rats  housed  In  Plexiglas  cages 
arrayed  on  the  tiers  of  a  stacked#  parallel-plate  circular  waveguide. 
Engineering  aspects  of  this  waveguide  and  the  exposure  environment  It  generated 
have  been  previously  reported  Cl],  The  sham-exposure  group  consisted  of  100 
cannulated  rats  housed  In  an  Identical,  but  unenergized#  collocated  facility. 
Results  reporting  blood  chemistry  and  hematology  In  these  same  animals  will  be 
published  In  the  next  volume  of  this  series.  Other  volumes  have  already 
published  results  on  adrenocorticotropic  hormone  (ACTH)  and  corticosterone  [2] 
and  prolactin  C3]. 

The  sympathetic-adrenal  medullary  system  plays  a  critical  role  In  the 
maintenance  of  cardiovascular  and  metabolic  homeostasis.  Plasma  catecholamines 
have  been  measured  to  assess  the  functional  activity  of  the  sympathetic-adrenal 
medullary  system  under  resting  conditions  or  during  stressful  stimulation. 
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Norepinephrine*  the  neurotransmitter  of  the  sympathetic  nervous  system* 
occurs  In  tissues  of  neural  crest  origin*  sympathetic  nerve  endings*  the  adrenal 
medulla*  and  other  chromaffin  tissues  as  well  as  In  the  brain.  Norepinephrine 
Is  synthesized  from  dopamine  by  the  enzyme  dopamlne-8 -hydroxylase  [41.  The 
predominant  sources  of  circulating  norepinephrine  are  sympathetic  nerve  endings 
and  the  adrenal  medulla.  Both  norepinephrine  and  dopamine- 6-hydroxylase  are 
secreted  from  sympathetic  nerve  terminals  In  proportional  amounts  during  nerve 
stimulation  [5*6]  and  can  be  accurately  measured  In  the  blood.  This  circulating 
norepinephrine  derives  largely  from  the  sympathetic  Innervation  to  vascular 
walls — especially  to  small  arteries  and  arterioles  which  provide  the  main  source 
for  peripheral  resistance  and  therefore  crucially  Influence  blood  pressure.  The 
extent  of  norepinephrine  "spillover"  from  the  synaptic  cleft  to  the  general 
circulation  depends  on  the  cleft  widths  perlsynaptlc  norepinephrine 
concentrations  are  relatively  low  for  narrow  gaps  but  high  for  wide  gaps  where 
the  concentrations  approach  those  estimated  to  be  attained  In  the  synapse. 
Since  vascular  Intramural  synapses  have  wide  gaps*  It  seems  likely  that  their 
proportional  contribution  to  circulating  norepinephrine  Is  large  when  compared 
to  nonvascular  noradrenergic  synapses  such  as  In  the  vas  deferens*  which 
typically  have  narrow  gaps.  Thus  the  level  of  plasma  norepinephrine  reflects 
both  adrenomedullary  and  sympathetic  nerve  activity. 

The  adrenal  medullary  responses  were  described  first  as  endosecretory 
responses  to  stress.  The  release  of  epinephrine  Is  part  of  this  response  and 
was  first  demonstrated  In  1914  [7]  In  cats  exposed  to  barking  dogs.  Similar 
responses  occur  after  many  psychological  or  physical  stimuli.  The  release  of 
epinephrine  correlates  with  the  degree  of  stress. 

The  physiologic  functions  of  the  dopamine  receptors  Include  vasodilation* 
Increased  sodium  excretion*  and  Increased  myocardial  contractility.  Even  change 
In  the  position  (from  standing  on  four  legs  to  exploring  the  cage  while  standing 
on  hind  legs)  Is  associated  with  enhanced  sympathetic  activity.  Similar  changes 
have  been  found  In  man  by  Sundln  [83. 

Exercise  Increases  plasma  catecholamines.  High  workloads  or  prolonged  work 
stimulates  several-fold  Increases  In  both  norepinephrine  and  epinephrine 
concentrations.  Many  other  stresses  Increase  the  release  of  catecholamines 
(particularly  epinephrine  and  norepinephrine).  Thus  the  plasma  level  of 
norepinephrine*  epinephrine*  and  (to  a  lesser  degree)  dopamine  fluctuates  widely 
In  a  mammal  reflecting  Increasing  or  decreasing  physical  activity  or  exposure  to 
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various  stressful  environments  [93.  The  determination  of  catecholamine  levels 
Is  used  to  quantitatively  measure  the  level  of  stress  Induced  on  the  autonomic 
nervous  system.  Sympathetic  neuronal  discharge#  with  adrenomedullary  release  of 
catecholamines  Into  the  blood#  Is  a  recognized  component  of  the  Immediate 
physiological  response  to  stress  [7»103.  Even  gentle  handling  produces  an 
Increase  In  epinephrine#  whereas  Immobilization  produces  massive  elevations  of 
circulating  levels  of  both  epinephrine  and  norepinephrine.  Decapitation  or 
restraint  lead  to  a  10-fold  Increase  In  circulating  norepinephrine  and  an  80- 
fold  Increase  In  circulating  levels  of  epinephrine#  whereas  dopamine  Increases 
to  a  lesser  degree  (Table  1).  The  high  levels  of  plasma  catecholamines  In  rats 
when  compared  with  other  animals  and  humans#  and  changes  produced  In 
pharmacological  and  physiological  experiments#  probably  reflect  environmentally 
Induced  changes  In  sympathoadrenomedullary  activity  rather  than  differences  In 
basal  sympathetic  neuronal  activity. 

TABLE  1.  CHANGES  IN  HORMONE  LEVELS  IN  CANNULATED  AND 
DECAPITATED  RATS 


Cannulated  After  decapitation 


Rat  # 

NOR 

EPI 

DA 

NOR 

EPI 

DA 

1 

825 

960 

185 

2 

104 

126 

30 

1275 

2795 

235 

3 

123 

144 

39 

1740 

1570 

210 

4 

144 

126 

25 

1870 

3565 

235 

5 

185 

113 

58 

1435 

2875 

170 

6 

174 

104 

76 

2660 

5430 

465 

7 

144 

159 

74 

1170 

1830 

365 

8 

153 

193 

28 

- 

- 

- 

9 

137 

154 

61 

1425 

2235 

205 

10 

162 

148 

74 

940 

1345 

260 

11 

144 

177 

43 

1520 

2975 

255 

12 

1930 

5295 

440 

X 

147 

144 

51 

1526 

2807 

275 

S.D. 

24 

28 

20 

515 

1485 

102 

All  hormone  concentrations  are  In  pg/mL. 
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II.  MATERIALS  AND  fCTHODS 


For  this  study#  the  concentrate ons  of  the  plasma  catecholamines 
norepinephrine#  epinephrine#  and  dopamine  were  chosen  as  sensitive  Indicators  of 
possible  environmental  stresses  Induced  by  RFR.  To  detect  and  quantitatively 
evaluate  possible  Increases  In  plasma  catecholamine  levels  induced  by  RFR#  blood 
was  sampled  and  assayed  from  65  exposed  and  64  sham-exposed  animals  (in  the  case 
of  epinephrine);  63  exposed  and  63  sham-exposed  animals  (In  the  case  of 
norepinephrine);  64  exposed  and  64  sham-exposed  animals  (In  the  case  of 
dopamine).  Analysis  of  the  data  obtained  from  the  blood  sample  assays 
determined  whether  there  were  any  RFR-Induced  changes  in  plasma  catecholamine 
concentrations. 

Animals.  The  rat  represents  a  comparatively  Inexpensive  and  homogeneous 
population.  For  this  reason#  It  Is  often  desirable  to  use  this  species  as  the 
animal  model  In  physiologic  studies. 

In  this  study#  male  Sprague-Dawley  rats  were  used.  All  experimental 
animals  were  obtained  from  the  same  building  and  room  at  CAMM  Research  Labs# 
Wayne#  New  Jersey.  The  animals,  weighing  approximately  60  g#  were  delivered  to 
Emory  University  where  they  were  caged  singly  and  given  water  and  food  (Purina 
Rat  Chow)  ad.  libitum.  Temperature  In  the  animal  rooms  was  maintained  at  24  +  1 
°C  and  the  photoperiod  was  12  hours/12  hours#  with  the  lighted  phase  occurring 
between  8  AM  and  8  PM. 

Experimental  Facility.  The  Georgia  Tech  Research  Institute's 
Radiofrequency  Radiation  Facility  [13  consisted  of  8  collocated  rooms  on  the 
basement  floor  of  the  Baker  Building  on  the  main  campus.  These  8  rooms  provided 
a  closed#  complete  facility  for  long-term  bioeffects  studies  Involving  rodents. 

The  100  exposure  and  100  sham-exposure  animals  were  housed  In  2  Identical# 
collocated  rooms  In  the  RFR  Facility.  Each  room  contained  a  stack  of  circular# 
parallel-plate  waveguides  fed  by  a  slotted-cyl Inder  antenna  system  for  radiating 
the  animals.  The  stacks  of  parallel  waveguides  consisted  of  five#  3.6-m  (12 
ft.)  diameter  plates  that  made  up  4  sets  of  circular  waveguides.  Twenty-five 
Individually  housed  rats  were  positioned  around  the  circumference  of  each 
waveguide  set.  The  walls  of  both  rooms  were  lined  with  anecholc  absorbing 
material  and  shielded  with  aluminum  foil  to  prevent  excessive  microwave  leakage 
radiation. 
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The  circular/  parallel-plate  waveguide  assembly  provided  a  1.0  mW/cnr 
exposure  field  around  the  circumference  of  the  plates.  The  45.7-cm  (18  In.) 
plate  separation  distance  permitted  propagation  of  a  TE^q  mode  wave  with 
horizontal  polarization.  The  power  density  displayed  a  cosine-squared 
dependency  between  the  plates,  with  the  maximum  power  density  occurring  midway 
between  each  set  of  plates.  This  arrangement  positioned  the  electric  field 
vector  parallel  to  the  rat’s  longitudinal  axis,  thereby  maximizing  the  coupling 
between  the  electric  field  and  the  rat. 

A  slotted-cyl  Inder  antenna  with  the  proper  diameter,  thickness,  slot 
length,  and  slot  width  dimensions  fed  the  stack  of  circular  waveguides  In  a 
manner  that  provided  an  essentially  constant  electric  field  Intensity  In  the 
azimuth  plane. 

Cages.  The  cages  were  constructed  of  clear  Plexiglas,  which  was 
essentially  RFR-transparent  at  435  MHz.  Clear  (rather  than  colored)  Plexiglas 
was  chosen  to  permit  visual  observation  of  the  rats.  Each  cage  was  22.9-cm  (9 
In.)  long  by  12.7-cm  (5  In.)  wide  by  17.8-cm  (7  In.)  tall.  These  dimensions 
complied  with  dimensions  recommended  by  the  National  Institutes  of  Health  for 
long-term  housing  of  rats  C1U.  The  food  hopper  and  water  bottle  were  placed  on 
the  distal  side  of  the  cage  to  minimize  their  Interaction  with  the  exposure 
field.  The  glass  floor  rods  In  the  cage  were  oriented  perpendicular  to  the 
cage’s  long  axis  to  Induce  the  rats  to  preferentially  align  themselves  parallel 
to  the  electric  field  vector.  The  slpper  tubes  of  the  water  bottles  were  made 
of  glass  to  be  nonperturbing  In  the  field.  Evaluations  of  the  cages  conducted 
In  the  circular,  parallel-plate  waveguide  assembly  showed  field  scattering  from 
the  Plexiglas  to  be  below  the  range  of  detection. 

The  RFR  Facility  contained  a  data  acquisition  system  for  storing  and 
processing  experimental  data,  an  electronic  balance  for  weighing  the  rats  during 
the  study,  and  rooms  for  transmitter  operation,  blood  sampling,  cage  washing, 
and  materials  storage. 

To  avoid  tfie  possible  effects  of  noise  during  this  study,  the  entire 
Radiation  Facility  was  kept  locked  to  avoid  unauthorized  entry.  Only  the  animal 
caretaker  and  the  technician  who  sampled  blood  from  the  animals  were  permitted 
uncontrolled  entry  to  the  Fatuity. 

Cannul atlon .  To  detect  and  quantitatively  evaluate  changes  In  plasma 
catecholamines,  the  resting  levels  of  these  hormones  first  had  to  be  determined. 
To  obtain  the  real  resting  values  of  the  three  hormones  In  undisturbed  animals. 


many  routine  techniques  for  handling  the  animals  and  for  sampling  the  blood  were 
unsuitable  for  this  study.  For  example,  guillotine  blood  sampling  techniques 
commonly  employed  in  many  endocrinological  studies  were  Immediately  ruled  out. 
To  use  each  animal  as  its  own  control,  arterial  blood  was  sampled  by  means  of 
chronically  Implanted  aortic  cannulas  [12,13,14].  This  simple,  inexpensive 
technique  permitted  remote,  stress-free  blood  sampling  in  conscious, 
unrestrained  and  resting  rats.  Arterial  blood  drawn  through  the  resting  rat's 
chronically  implanted  cannula  was  assayed  for  plasma  norepinephrine,  plasma 
epinephrine,  and  plasma  dopamine. 

The  idea  of  sampling  venous  blood  from  the  animals  was  abandoned.  In 
venous  blood  vessels,  the  flow  regime  is  laminar  with  blood  flowing  In  discrete 
layers.  The  layers  of  blood  in  the  middle  of  the  vessels  travel  much  faster 
than  those  close  to  the  vessel  walls.  The  most  Important  consideration, 
however,  was  that  blood  layers  do  not  mix  in  venous  blood  vessels.  Thus,  a 
sample  of  venous  blood,  withdrawn  with  a  needle  or  a  cannula,  might  represent 
the  blood  returning  from  one  part  of  the  body  or  the  other,  from  a  single  organ 
or  muscle,  or  from  any  one  of  the  endocrine  glands.  For  this  reason,  we  decided 
to  sample  arterial  blood,  which  Is  always  fully  mixed.  The  mixing  occurs  In  the 
left  ventricle  of  the  heart  and  In  early  parts  of  the  aorta.  Only  small  amounts 
of  arterial'  blood  (up  to  0.6  mL)  were  withdrawn  from  resting  rats  about  once 
every  3  to  5  weeks.  Removing  greater  volumes  of  blood  has  been  shown  to  elevate 
plasma  norepinephrine  concentrations  In  the  rats  (Fig.  1). 

We  used  PE-10  arterial  cannulas  In  this  study.  Larger  PE-50  cannulas  were 
unsuitable  because  they  could  develop  large  blood  clots  if  not  drained 
frequently.  Large  cannulas  require  multiple  flushing  to  remain  patent,  but 
flushing  might  Induce  multiple  strokes  In  the  animals.  Chronic  cannulation  of 
the  aorta  with  a  PE-10  cannula  was  preferable  to  cannulation  of  other  arterial 
blood  vessels.  Cannulation  of  the  abdominal  aorta  provided  long-term  functional 
cannulas,  but  the  cannulation  procedure  was  lengthy  (20-30  min)  and  required 
opening  the  abdominal  cavity  and  temporary  dislocation  of  the  gastro-intestinal 
system.  The  abdominal  aortic  cannula  had  a  much  larger  dead  space  than  the 
aortic  cannula.  Cannulation  of  the  aorta  through  the  left  carotid  artery,  on 
the  other  hand,  required  an  incision  of  1-1.5  cm  that  neither  penetrated  body 
walls  nor  entered  the  abdominal  cavity.  Further,  this  cannulation  could  be 
completed  in  about  8  min. 
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Figure  1.  Effect  of  1.0  mL  bleeding  on  resting  plasma  norepinephrine 
concentration. 


The  carotid  artery  of  the  animal  was  cannulated  8  to  10  days  before  the* 
animals  entered  the  study.  The  surgery  was  done  using  ketamlne-xylazlne 
anesthesia  (1:1  mixture;  ketamine  100  mg/mL,  xylazlne  20  mg/mL,  l.m.  0.1  mL  / 
100  g  of  body  weight).  The  catheter  was  filled  with  slightly  heparinized 
saline*,  and  the  distal  end  was  sealed  with  a  nylon  plug.  Stress  hormone  levels 
returned  to  the  basal  values  about  3  days  after  Implantation  of  the  chronic 
arterial  cannulas.  The  first  blood  sampling  occurred  10  days  after  aortic 
cannulatlon. 

Blood  Sampling.  Although  the  half-life  of  plasma  catecholamines  Is  only  1 
to  3  min  C15],  a  strong  stimulus  leaves  plasma  catecholamine  levels  relatively 
high  for  a  period  of  up  to  15-20  min.  Normal  handling  (lifting  the  rat)  evoked 
a  75%  Increase  In  epinephrine  concentration  accompanied  by  a  small  Increase  in 
norepinephrine  concentration.  However,  the  animals  had  to  be  handled  when  they 
were  removed  from  their  exposure  cage  and  placed  In  the  "sampling  box"  In 
preparation  for  blood  withdrawal.  To  avoid  the  undesired  effects  of  handling  on 
catecholamine  levels,  blood  from  the  aortic  cannula  was  sampled  30  min  after  the 
animal  was  placed  In  the  sampling  box.  This  procedure  permitted  the  altered 
plasma  catecholamine  levels  sufficient  time  to  return  to  their  basal  (resting) 


^,5-cm3 heparin  sodium  (from  beef  lung),  1000  unlts/mL  per  30  cm3 saline. 


values.  Each  animal  was  preconditioned  for  the  sampling  box  through  a  regime  of 
several  30-mln-long  experiments  conducted  during  a  1-week  period  before  entering 
the  study. 

After  acclimating  for  30  min  In  the  sampling  box#  the  rat’s  cannula  was 
positioned  through  the  slot  In  the  top  of  the  box  (Fig.  2).  The  heparinized 
saline  was  then  removed  from  the  cannula#  and  a  0.6  mL  blood  sample  was  taken 
from  the  resting  rat  using  a  sterile  1-cm3 tubercul In  syringe  fitted  with  a  30- 
ga  needle.  The  syringe  and  needle  were  rinsed  with  ethylene  glycol-bls 
tetraacetlc  acid  (EGTA)/glutath1one  before  sampling.  The  blood  sample  was 
placed  In  an  EGTA/gl utathlone-treated  1.0  mL  capillary  blood  collection 
container  (prepared  In-house  and  stored  under  refrigeration  to  prevent  chemical 
breakdown)#  shaken#  and  then  placed  on  Ice.  The  blood  sampling  procedure 
required  about  2  min  for  each  rat. 


Figure  2.  Sampling  blood  from  the  chronic  aortic  cannula  of  a  resting# 
unrestrained#  and  unanesthetized  rat. 

Plasma  catecholamine  levels  In  conscious  unrestrained  rats  with  chronic 
Indwelling  catheters  were  considerably  lower  than  previously  reported  for  the 
rat  [16]. 


31coa  Sampling  Scheeule.  Figure  3  shows  the  sampling  schedule  designee  for 
the  experiment.  The  200  rats  were  introduced  into  the  study  in  4  groups  of  SO 
animals  each.  The  groups  entered  in  a  staggered  manner  to  facilitate  the 
process  of  logging-in  and  establishing  the  new  animals.  Each  group  contained  25 
exposure  ana  25  sham-exposure  animals.  Of  the  25  exposure  (or  sham-expesure) 
animals,  20  were  sampled  for  plasma  stress  hormones,  while  the  remaining  5  were 
used  for  hematology  studies. 

The  sampling  duration  was  36  weeks  long.  Including  a  6-week  preexposure 

adaptation  period,  a  24-week  exposure  period,  and  a  6-week  pcstexpcsure  period. 

With  group  staggering  taken  into  account,  the  experiment  duration  was  42  weeks 

long  (since  introduction  of  the  4  groups  was  staggered  in  2-week  ■intervals). 

♦ 

Plasma  catecholamines  were  to  be  sampled  for  all  periods  marked  (B)  in  Figure  3. 
Therefore,  each  animal  should  have  been  sampled  for  plasma  norepinephrine, 
epinephrine,  and  dopamine  at  weeks  -5,  -2,  1,  4,  7,  ...,  28.  This  schedule  was 
rather  rigorous  and  therefore  could  tolerate  slight  fluctuations  in  protocol 
without  ill  effects. 
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Figure  3.  Sampling  and  exposure  timetable. 


Plasma  Catecholamine  Determinations.  Plasma  catecholamines  were  measured 
with  a  radloenzymatlc  method  according  to  Penler  and  Johnson  [17].  Briefly#  the 
three  catecholamines  were  first  converted  to  their  o-methylated  analogues  by 
catechol-o-methyl-transferase  In  the  presence  of  S-adenosyl-methomlne-^H  and 
thereafter  extracted  following  addition  of  sodium  tetraphenyl byrate.  This 
extraction#  together  with  an  improved  quick  chromatographic  separation  and  the 
oxidation  of  the  epinephrine  and  norepinephrine  derivatives. to  vanillin,  yielded 
an  extremely  high  sensitivity  and  specificity  of  the  method.  The  assay  allowed 
the  determination  of  norepinephrine#  epinephrine,  and  dopamine  In  plasma  volumes 


of  20-100  yl. 


III.  RESULTS  AND  ANALYSIS 


Plasma  Norepinephrine.  Appendix  A  contains  the  data  collected  during  the 
preradiation  and  radiation  periods  for  both  the  exposure  and  sham-exposure 
groups.  The  high  variance  displayed  by  the  data  for  the  entire  sampling  period 
indicated  various  degrees  of  animal  activity  at  the  time  of  blood  sampling. 
Since  the  boxes  had  opaque  walls#  the  activity  of  each  animal  before  sampling 
was  not  recorded.  However#  as  previously  mentioned#  it  was  unlikely  that  the 
stimulation  of  placing  the  rats  in  the  sampling  boxes  had  a  major  effect  on 
resting  norepinephrine  concentration#  since  the  increase  in  norepinephrine 
secretion  induced  by  animal  handling  would  disappear  20  to  30  min  following  the 
stress. 

Figures  4  and  5  present  the  raw  norepinephrine  concentration  in  scatter 
diagram  form  (the  dotted  lines  pass  through  the  mean  response  at  each  week  data 
were  collected).  Despite  a  3 -week  effort  to  precondition  the  animals  to  the 
sampling  box  environment  before  drawing  blood  samples#  the  basal  resting  value 
of  plasma  norepinephrine  decreased  during  weeks  -3#  -2#  and  0.  This  same 
behavior  was  also  observed  in  plasma  A CTH#  plasma  corticosterone,  and  plasma 
prolactin  [2,3].  After  the  first  week#  the  data  displayed  a  nearly  linear 
response.  The  "spikes"  occurring  at  weeks  10  and  17  (sham-exposure  group)  are 
the  mean  values  resulting  from  7  and  3  observations#  respectively;  the  spike  at 
week  11  (exposure  group)  is  the  mean  value  resulting  from  5  observations.  The 
wide  range  spanned  by  the  2-slded  95/5  confidence  Interval  at  each  value 
indicated  that  these  "spikes"  may  not  represent  drastic  deviations  from  the 
established  norepinephrine  resting  concentration.  Noise  and  unfamiliar  persons 
visiting  the  Radiation  Facility  may  have  also  contributed  to  the  sham-exposure 
group  spike  at  week  10. 

Mean  plasma  norepinephrine  concentrations  in  the  exposure  and  sham-exposure 
groups  did  not  appear  significantly  different  when  plotted  on  the  same  axis 
(Fig.  6).  This  was  preliminary  evidence  Indicating  that  chronic  exposure  to 
435-MHz  RFR  did  not  affect  the  resting  level  of  plasma  norepinephrine.  A 
statistical  analysis  was  subsequently  performed  on  the  data  to  test  this 
hypothesis. 

The  analysis  Involved  using  multiple  linear  regression  techniques  to  build 
a  model  describing  plasma  norepinephrine  levels  as  a  function  of  time  and 
incident  RFR.  'arms  of  the  polynomial  model  thus  obtained  were  tested  for  their 
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Figure  4.  Norepinephrine  concentration  data  scatter  diagram  (sham-exposure  group). 
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Figure  5.  Norepinephrine  concentration  data  scatter  diagram  (exposure  group). 
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Figure  6.  Mean  plasma  norepinephrine  concentrations  versus  time. 


significance  in  describing  the  collected  data.  Various  diagnostic  procedures# 
including  model  lack-of-fit  tests#  residual  analysis#  and  autoregressive 
analysis#  were  then  applied  to  the  model  to  check  Its  validity.  Appendix  B 
contains  a  detailed  description  of  this  statistical  methodology#  as  well  as  the 
individual  analyses  for  each  of  the  three  catecholamines. 

The  statistical  analysis  Indicated  that  there  was  no  significant  difference 
between  the  sham-exposure  and  exposure  groups.  The  final  polynomial  model  was 
solely  a  function  of  time.  Resting  norepinephrine  levels  were  at  their  highest 
value  (approximately  299  pg/mL#  as  calculated  from  the  model  derived  In  the 
norepinephrine  statistical  analysis  of  Appendix  B)  at  the  study  onset  (week  -3). 
The  resting  level  then  gradually  declined#  reaching  Its  lowest  point  of  an 
estimated  222  pg/mL  at  week  13  of  the  study.  Norepinephrine  concentration  then 
appeared  to  rise#  reaching  a  value  of  about  232  pg/mL  at  week  19  of  the  study# 
which  was  the  last  week  for  which  data  was  available.  Since  no  data  were  taken 
beyond  week  19,  there  was  no  effort  to  extrapolate  a  value  for  week  29  of  the 
study. 

Further  analysis  determined  the-  smallest  change  In  resting  norepinephrine 
concentration  (between  exposure  and  sham-exposure  groups)  that  the  protocol  was 
capable  of  detecting.  If  there  were  any  RFR-Induced  effects  on  the  resting 
concentration  of  norepinephrine#  they  would  have  to  lie  within  the  range  of  +  15 
pg/mL  from  the  estimated  resting  concentration  of  273  pg/mL.  Since  values  of 
norepinephrine  between  258  pg/ml  and  288  pg/mL  are  considered  normal  In 
unstressed  rats#  there  was  no  Indication  that  chronic  RFR  exposure  resulted  In 
any  stress  to  the  animals#  as  measured  by  plasma  norepinephrine. 

Plasma  epinephrine.  Appendix  G  contains  the  data  collected  during  the  pre¬ 
radiation  and  radiation  periods  for  both  exposure  and  sham-exposure  groups. 
Like  norepinephrine#  this  hormone  also  displayed  a  variance  about  the 
established  resting  level  due  to  varying  amounts  of  animal  activity.  Since 
plasma  epinephrine  concentrations  were  sensitive  to  handling  and  related 
stresses#  each  animal  was  given  30  min  to  allow  the  epinephrine  concentration  to 
return  to  the  basal  value. 

Figures  7  and  8  present  the  raw  epinephrine  concentration  data  in  scatter 
diagram  form  (the  dotted  lines  pass  through  the  mean  epinephrine  response  at 
each  week  data  were  collected).  Once  again#  the  mean  epinephrine  values  In  both 
exposure  and  sham-exposure  groups  declined  In  the  Initial  3  weeks  of  the  study. 
This  decline  was  attributed  to  the  animals  being  Inadequately  preconditioned  to 
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the  sampling  boxes.  Once  the  animals  adapted  to  the  sampling  box  environment* 
the  epinephrine  concentrations  In  both  RFR-exposed  and  sham-exposed  animals 
remained  about  the  same.  The  small  amount  of  "splklness"  In  the  plots  was  the 
random  effect  of  sampling  within  a  population. 

The  mean  epinephrine  concentrations  in  the  exposure  and  sham-exposure 
groups  did  not  seem  to  be  significantly  different  when  the  two  plots  were 
compared  to  one  another  (Fig.  9).  This  evidence  suggested  tnat  chronic  exposure 
to  435-MHz  RFR  did  not  affect  the  resting  concentrations  of  plasma  epinephrine. 
A  statistical  analysis  was  then  performed  on  the  epinephrine  data  to  test  this 
hypothesis. 

The  statistical  analysis  Involved  building  a  polynomial  function  relating 
epinephrine  concentration*  time,  and  RFR  radiation  In  the  same  manner  as  the 
previous  hormones  (ACTH,  corticosterone,  prolactin,  and  norepinephrine).  The 
terms  of  the  polynomial  model  were  then  tested  to  determine  their  significance 
in  describing  the  epinephrine  data  set.  The  final  model,  which  was  Independent 
of  RFR,  was  then  verified  using  1  ack-of-f  It,  residual  analysis*  and 
autoregression  techniques.  The  complete  statistical  analysis  Is  Included  In 
Appendix  B. 

The  analysis  concluded  that  RFR  had  no  effect  on  the  exposure  group  when 
compared  to  the  sham-exposure  group.  Epinephrine  concentration  during  the  study 
did  display  a  time  dependence,  however, '  decreasing  from  an  estimated  Initial 
concentration  of  181  pg/mL  at  the  study  onset  (week  -3)  to  a  low  of  119  pg/mL 
during  the  exposures  (week  12),  and  then  Increasing  to  about  134  pg/mL  at  week 
19,  the  last  week  for  which  data  were  available.  Once  again,  no  effort  was  made 
to  use  the  epinephrine  model  as  a  forecasting  tool  for  week  29.  Further 
analysis  Indicated  that.  If  there  were  any  RFR-induced  effects,  they  had  to  lie 
within  a  range  of  +.13  pg/mL  from  the  resting  value  of  159  pg/mL.  Since  resting 
epinephrine  concentrations  between  146  pg/ml  and  172  pg/mL  are  considered  normal 
In  unstressed  rats,  there  was  no  Indication  that  long-term  RFR  exposure  produced 
any  stress  as  measured  by  plasma  epinephrine  concentrations. 

Plasma  dopamine.  Appendix  L  contains  the  data  collected  during  the  pre¬ 
radiation  and  radiation  periods  for  both  exposure  and  sham-exposure  animals. 
The  variance  In  the  data,  as  mentioned  before,  derived  principally  from  various 
levels  of  animal  activity  Immediately  before  sampling.  The  30-mln  acclimation 
time  allowed  dopamine  concentrations  to  return  to  the  resting,  basal  level. 


Figures  10  and  11  present  the  raw  dopamine  concentration  data  in  scatter 
diagram  form  (the  dotted  lines  pass  through  the  mean  dopamine  response  at  each 
week  data  were  collected).  Again#  the  mean  aopamine  response  in  both  exposure 
and  snam-exposure  groups  declined  in  the  initial  3  weeks  of  the  study.  This 
decline  was  similar  to  the  observations  noted  in  the  other  hormones  assayed 
(ACTW,  corticosterone#  prolactin#  norepinephrine#  and  epinephrine)  anc  w  as 
attributed  to  the  same  source  (animals  inadequately  preconditioned  to  sampling 
box).  After  the  animals  adapted  to  the  sampling  box  environment#  the  dopamine 
values  in  both  exposure  and  sham-expcsure  groups  tended  to  stabilize  (weeks  1 
through  19).  The  spikiness  in  the  plots  was  a  result  of  random  sampling  within 
both  populations. 

Mean  plasma  dopamine  concentrations  did  not  appear  to  be  larger  in  the 
exposure  group  when  compared  to  the  sham-exposure  group  (Fig.  12).  If  anything, 
the  opposite  seemed  to  be  the  case  for  the  length  of  the  experiment.  _  This 
result  indicated  that  chronic  exposure  to  435-MHz  RFR  did  not  induce 
physiological  changes  in  the  rat  population  that  were  manifested  as  increased 
resting  dopamine  concentrations.  A  statistical  analysis  was  therefore 
performed  on  the  data  set  to  test  this  hypothesis. 

The  statistical  analysis  performed  was  identical  in  procedure  to  that  used 
in  the  analysis  of  the  other  study  hormones.  A  detailed  description  of  the 
general  methodology#  and  the  specific  dopamine  analysis#  is  given  in  Appendix  B. 
The  analysis  for  all  hormones  used  the  SAS  Statistical  Software  resident  on  the 
Georgia  Tech  IBM  4381  mainframe  to  run  tests  and  produce  the  analysis  hardcopy. 

The  analysis  gave  no  indication  of  Increased  plasma  dopamine  in  the 
exposure  group  when  compared  to  the  sham-exposure  group.  In  fact#  the  estimated 
dopamine  concentration  in  the  exposure  group  remained  significantly  smaller  than 
that  of  the  sham-exposu re  group  from  the  initiation  of  exposures  to  the 
termination  of  the  experiment.  Resting  dopamine  values  were  at  their  highest 
for  week  -3  of  the  study  (about  62  pg/mL  sham-exposed#  65  pg/mL  exposed).  The 
resting  levels  of  both  groups  then  declined#  reaching  the  lowest  value  of  32 
pg/mL  at  week  12  (sham-exposed);  20  pg/ml  at  week  16  (exposed).  Beyond  these 
points#  dopamine  concent  rat  1  on  gradually  Increased#  with  estimated 
concentrations  of  39  pg/ml  (sham-exposed)  and  21  pg/mL  (exposea)  at  week  19#  the 
last  week  for  which  data  were  collected. 

Further  analysis  showed  that  the  smallest  change  in  resting  dopamine 
concentration  that  the  protocol  could  reliably  detect  was  about  6  pg/ml  above  or 


below  an  estimated  resting  concentration  of  51  pg/mL.  This  analysis  indicated 
that,  if  RFR  increased  resting  dopamine  levels  above  57  pg/ml,  the  protocol 
would  have  found  a  significant  positive  RFR  effect.  In  fact,  dopamine 
concentrations  of  up  to  120  pg/mL  were  considered  normal-  for  a  population  of 
healthy,  unstressed  Sprague-Dawley  rats.  Therefore,  there  was  no  indication 
that  chronic  exposure  to  low-level  435-MHz  RFR  produced  any  stress  in  the 
exposure  group  (when  compared  to  the  sham-exposure  group)  as  measured  by  ‘the 
concentration  of  blood-borne  dopamine. 
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IV.  DISCUSSION 


Minute  amounts  of  free  (unconjugated)  catecholamines  are  normally  found  In 
both  human  and  animal  blood  plasma.  These  hormones  undergo  rapid  changes  which 
reflect  sympathetic  nerve  activity  [18,19].  The  radloenzymatic  techniques 
available  for  quantitative  determinations  of  norepinephrine,  epinephrine,  and 
dopamine  in  a  few  microliters  of  plasma  permit  monitoring  of  the  sympathoadrenal 
activity  In  small  laboratory  animals  such  as  the  rat. 

Arterial  blood  pressure,  ambient  temperature,  body  temperature, 
physiological  activity,  and  certain  biological  characteristics  (e.g.»  animal 
strain)  have  an  effect  on  the  level  of  circulating  plasma  catecholamine 
concentrations  [20],  Different  strains  of  rats  have  dissimilar  levels  of 
resting  catecholamines  [21,22].  Both  normotenslve  and  hypertensive  rats  show 
the  same  catecholamine  response  at  rest,  but  hypertensive  rats  show  a  greater 
catecholamine  response  during  stress  [21,23]. 

For  a  particular  strain  of  animal,  the  resting  level  of  plasma 
catecholamines  Is  always  the  same  [24],  permitting  measurement  of  increases  in 
plasma  catecholamine  concentration  and  (from  these  Increases)  evaluation  of  the 
level  of  stress  an  animal  underwent  [25,26].  The  stronger  a  stress  and  the 
longer  Its  duration,  the  higher  the  concentration  of  plasma  epinephrine, 
norepinephrine,  and  In  some  cases  dopamine  [27],  Even  a  small  ■  reduction  of 
blood  volume  Increases  plasma  catecholamine  levels  [28,29]. 

To  obtain  reliable  measurements  of  circulating  catecholamines  In  rats 
required  appropriate  methods  for  blood  collection  to  avoid  catecholamine 
Increase  due  to  physical  stress  [16],  In  this  study,  resting  levels  of 
catecholamines  were  considerably  lower  In  rats  whose  blood  samples  were 
collected  from  Indwelling  cannulas  than  values  where  blood  was  obtained  by 
decapitation  or  other  stressful  methods. 

The  results  of  our  experiment  Indicate  that  exposure  to  chronic  low-level 
RFR  did  not  represent  a  stress  measurable  as  an  Increase  In  norepinephrine  and 
epinephrine  concentration  of  Irradiated  rats.  Similar  results  were  obtained 
when  plasma  ACTH,  plasma  corticosterone,  and  plasma  prolactin  were  determined  In 
Identical  situations  [2,3]. 

In  this  study,  plasma  dopamine  decreased  In  RFR-exposed  animals.  Though 
significant,  the  small  plasma  dopamine  decrease  might  not  be  physiologically 
Important.  It  would  be  of  Interest  to  ascertain  whether  this  lowered  dopamine 
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concentration  persists  after  RFR  exposure  is  Interrupted  for  several  cays  or 
weeks  (the  rats  were  removed  from  the  RFR  field  for  30  min  to  obtain  the  blood 
samples).  The  large  individual  variation  observed  in  the  plasma  catecholamine 
levels  of  both  RFR-exposed  and  sham-exposed  animals  was  probably  the  consequence 
of  various  levels  of  animal  physiological  activity  during  or  just  before  blood 
sampling.  It  is  known,  for  instance,  that  during  sleep  plasma  levels  of 
norepinephrine  and  epinephrine  are  below  those  of  the  resting  state  [30]. 

Although  plasma  catecholamine  half-life  is  only  1  to  3  min  [15],  a  strong 
stimulus  leaves  plasma  catecholamine  levels  relatively  high  for  10  to  15  min. 
For  this  reason,  blood  was  sampled  from  the  resting  animals  30  min  after  gentle 
placement  in  tne  sampling  boxes,  permitting  plasma  catecholamine  levels  to 
return  to  the  resting  level. 

During  the  6-month  study  duration,  the  rats  aged  somewhat.  Some 
investigators  have  reported  changes  In  catecholamine  secretion  induced  by  aging 
[31,32,33].  However,  new  studies  demonstrate  that  aging  does  not  change  the 
rat’s  responsiveness  to  either  internal  or  external  stimuli  that  evoke 
catecholamine  secretion  [34].  The  same  study  failed  to  find  changes  over  a 
several  month  period  in  resting  plasma  catecholamine  concentration  of  rats. 

In  conclusion,  our  results  indicated  that  a  435-MHz  pulsed-wave  environment 
did  not  Increase  resting  plasma  catecholamine  concentrations  In  rats.  The 
statistical  analysis  of  the  data  Indicated  that  If  there  were  any  RFR-Induced 
effects  on  resting  plasma  catecholamine  concentrations,  they  would  lay  within  a 
range  of  +  15  pg/mL  from  an  estimated  resting  concentration  of  273  pg/mL  in 
norepinephrine;  +.  13  pg/ml  from  an  estimated  resting  concentration  of  159  pg/mL 
in  epinephrine;  and  +.  6  pg/mL  from  an  estimated  resting  concentration  of  51 
pg/mL  In  dopamine.  These  values  are  not  typical  of  rats  exposed  to  stress. 
Therefore,  this  study  concludes  that  a  1.0  mW/cm^  435-MHz  pulsed-wave  (1.0  ys 
pulse  width,  1  kHz  pulse  rate)  RFR  environment  did  not  induce  any  detectable 
increase  in  stress,  as  measured  by  resting  catecholamine  concentrations  in  the 
exposure  group  of  cannulated  male  Sprague-Dawley  rats  when  compared  to  the  sham- 
exposure  group. 
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STATISTICAL  METHODOLOGY 

The  balanced  design  of  this  experiment  (requiring  that  25  animals  from  each 
100  animal  group  be  sampled  once  every  3  weeks  for  stress  hormones)  should  have 
produced  data  easily  tested  by  balanced*  2-way  analysis  of  variance  (ANOVA) 
statistics  with  12  levels  of  factor  A  (time)  and  2  levels  of  factor  3  (RF 
raciation).  However*  data  collection  did  not  proceed  according  to  protocol  in 
that*  in  numerous  cases*  samples  were  collected  at  odd  intervals  ( inval i cati no 
the  orthogonality  of  the  design)  and  the  number  of  samples  "taken  per  week  varied 
above  and  below  the  25  animal  mark  (unbalancing  the  design).  These  two  factors 
combined  to  lower  the  power  of  ANOVA  statistics  (power  being  defined  as  the 
ability  to  reject  the  null  hypothesis  given  the  null  hypothesis  should  be 
rejected)  trying  to  test  the  model 

y)jk  =  K  T,  +  Bj  +  refJ  +  «,jk.  (B-l) 

hormone  concentration  (response)* 

the  normal  hormone  resting  concentration, 

the  change  in  hormone  resting  concentration  induced  by  RFR, 

the  change  in  hormone  resting  concentration  induced  by  time, 

the  change  in  hormone  resting  concentration  induced  by  the 
Interaction  between  RFR  and  time,  and 

noise  within  the  system  (sampling  and  assaying  errors) 
for  the  following  hypotheses: 

H0:  T0  =  T1  =  °’ 

Hj:  Tq  or  Tj.  A  0  (RFR-induced  effects),  (B-2) 

H0:  =  62  =  .  .  .  =  3. .  =  0. 

at  least  one  Sj  4  Q  ftime-induced  effects),  (B-3) 

Hg :  tS  =  0*  and 

H^:  at  least  one  4  0  (interaction  between  RFR  and  time).  (3-4) 


where 


y 


i  jk 
V 


Ti = 

3,.  = 


tS 


ij 


'1  jk 
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However#  examination  of  the  collected  data  suggested  an  alternative 
approach  in  that  the  data  resembled  what  might  have  been  collected  in  an 
unplanned  experiment  monitoring  over  time  the  operation  (in  this  case# 
characterized  by  resting  animal  hormone  concentrations)  of  an  established  RF 
radiation  facility.  Data  of  this  type  are  often  successfully  treated  by 
employing  linear  regression  techniques  to  develop#  build#  and  test  a  linear  (or 
intrinsically  linear)  model  whose  parameters  can  be  used  to  predict  the  system 
response  at  various  treatment  levels.  Therefore#  we  decided  to  proceed  with  a 
regression  approach  to  data  analysis. 

Plasma  Norepinephrine  Statistical  Analysis. 

Examination  of  the  norepinephrine  scatter  diagrams  of  Figures  4  and  5  yield 
an  essentially  linear  norepinephrine  response  versus  time  beyond  week  0  of  the 
study.  There  was,  however,  a  certain  amount  of  positive  curvature  present  at 
both  the  study  initiation  and  study  conclusion#  particularly  in  the  sham- 
exposure  group.  Therefore#  a  quadratic  polynomial  function  was  empirically 
chosen  to  test  for  RFR  effects  within  the  exposure  and  sham-exposure  groups. 
Thus#  the  norepinephrine  response  was  modelled  with  a  nonzero  Intercept  6g  and 
an  RFR-induced  effect  on  this  Intercept  (ogz),  a  nonzero  linear  slope  8^  and  an 
RFR-Induced  effect  on  this  slope  (^z)#  and  a  quadratic  coefficient  3  and  RFR- 
induced  effect  on  this  curvature  (ct^z).  The  statistical  significance  of  these 
terms  determined  the  Importance  of  their  contribution  to  the  final  model.  The 
equation  describing  the  initial  model  was  therefore: 

y  =  s0  +  +  3gZ  +  ^zx  +  anzx^  (B-5) 

where  y  =  plasma  norepinephrine  concentration  (In  pg/mL)# 
x  3  time  (In  weeks)#  and 

z  3  a  categorical  variable  with  value  0  for  animals  in  sham-exposure 
group  and  value  1  for  animals  In  exposure  group. 

Raw  data  from  the  norepinephrine  spreadsheet  (Appendix  A)  were  put  on 
computer  file.  A  Statistical  Analysis  System  (SAS)  formatting  program  (Appendix 
C)  was  prepared  to  read  the  data  and  perform  the  desired  statistical  tests  on 
the  model . 

The  first  test  Identified  terms  within  the  model  which  contributed  the 
least  toward  forming  a  statistically  significant  regression.  These  procedures 
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were  used  in  comDination  with  an  Initial  regression  on  the  general  model  (not 
included)  to  evaluate  the  statistical  significance  of  terms  modelling  the 
norepinephrine  concentration  time  dependency  and  terms  modelling  the  RFR-induced 
effects  on  norepinephrine  concentration.  Two  types  of  model  "building” 
procedures  were  used:  forward  stepwise  regression  and  maximum  R^  regression. 
Forward  stepwise  regression  produced  a  model  by  calculating  F  statistics  for  all 
variables  not 'in  the  model  #  ar.d  then  adding  a  variable  to  the  model  if  its  F 
statistic  was  significant  at  a  given  a  risk  (for  this  reason#  the  forward 
procedure  begins  with  no  variables  in  the  model).  Once  a  variable  was  added  to 
the  model#  the  procedure  recalculated  F  statistics  for  all  the  terms  in  che 
model#  and  rejected  any  terms  whose  F  statistic  rose  above  a  given  a  risk.  In 
this  manner#  forward  stepwise  regression  eventually  settled  on  a  model  including 
all  terms  whose  3  risk  was  low  enough  to  permit  initial  entry  and  then  not  be 
rejected  upon  the  addition  of  other  terms. 

Maximum  R^  regression  took  this  procedure  further,  producing  lists  of  the 
best  1-parameter  model#  best  2-parameter  model,  best  3-parameter,  etc.,  until 
all  of  the  parameters  were  Included  in  the  final  model.  This  procedure 
permitted  discrimination  of  different  models  using  number  of  parameters  as  a 
judgement  criterion. 

Both  forward  stepwise  and  maximum  R  regressions  indicated  that  the  model 
which  best  fit  the  data  was: 


where 


y  = 


11 


x2# 


s0  =  272.8, 


8ii " 


-7.79, 

0.30. 


and 


(5-6) 


The  entry  and  exit  a  risk  was  0.10.  The  outputs  of  both  regression 
procedures  are  Included  In  Appendix  D.  Note  that  the  absence  of  2  terms 
indicated  that,  at  a  0.10  risk#  there  was  no  statistical  difference  in  plasma 
norepinephrine  concentrations  between  the  exposure  and  sham-exposure  group.  The 
estimated  resting  concentration  of  plasma  norepinephrine,  272.8  pg/mL  (3q), 
agreed. well  with  established  values  cited  in  the  literature  (300  +40  pg/mL). 
This  agreement  was  an  indication  of  no  systematic  error  within  the 
sampling/assaying  procedure. 
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3oth  exposure  and  sham-exposure  groups  did  display  a  time  dependency  in 
norepinephrine  concentration.  Resting  norepinephrine  levels  were  at  their 
highest  value  (about  298.9  pc/mL)  at  the  study  onset  (week  -3).  The  resting 
level  then  gradually  declined*  reaching  its  lowest  point  of  221.8  pg/mL  at  week 
13  of  the  study.  Norepinephrine  concentration  then  seemed  to  rise,  reaching  a 
value  of  232.0  pg/ml  at  week  19  of  the  study,  which  was  the  last  week  data  were 
taxen.  All  of  the  just  mentioned  values  were  well  within  the  normal  bounds  of 
plasma  norepinephrine  concentration  in  healthy,  unstressed  rats.  Therefore,  it 
seemec  that  chronic  exposure  to  435-MHz  RFR  did  not  result  in  an  increase  in 
stress  (as  measured  by  the  concentration  of  plasma  norepinephrine)  in  the 
exposure  group  when  compared  to  the  sham-exposure  group. 

Tne  just  mentioned  conclusions  could  only  be  accepted  once  the  assumptions 
used  to  Duild  the  final  model  were  verified.  These  assumptions  included  nc 
model  lack-of-fit,  NID(0»  a^)  resiaual  distribution  (meaning  residuals  were 
normal  anc  independently  distributed  with  mean  zero  and  variance  o2),  and  no 
mocel  mult icoll inear ity. 

Since  multiple  observations  of  norepinephrine  concentration  were  taken  for 
tne  weeks  containing  data,  it  was  possible  to  perform  a  model  lack-of-fit  test 
on  the  regression.  The  lack-of-fit  involved  breaking  the  sum-cf-squares  errcr 
from  the  regression  into  two  components:  sum-of-squares  pure  error, 
representing  the  actual  variation  due  to  the  sampling  and  assaying  process  and 
sum-of-squares  lack-of-fit,  representing  the  variation  due  to  the  difference 
between  the  mean  value  at  one  week  when  compared  to  the  fitted  value  at  the  same 
week.  A  test  statistic  was  then  computed  comparing  the  sum-of-squares  lack-of- 
fit  to  the  sum-of-squares  pure  error;  sufficiently  high  values  of  the  test 
statistic  indicated  model  lack-of-fit. 

Sum-of-squares  error  was  obtained  from  the  ANCVA  table  produced  in  the 
regression  procedure  output.  Sum-of-squares  pure  error  was  obtained  by 
analyzing  the  experiment  from  2-way,  fixed  effects  ANOVA  viewpoint.  The  sum-of- 
squares  lack-of-fit  was  then  computed  from  the  difference  of  sum-of-squares 
error  minus  sum-of-squares  pure  error.  Calculations  to  compute  the  critical 
value  Fg  are  detailed  in  Appendix  E. 

Since  the  computed  test  statistic  Fg  was  smaller  than  the  critical  value, 
there  was  insignificant  model  lack-of-fit.  This  indicated  that  the  quadratic 
function  modelling  norepinephrine  concentration  versus  time  was  a  good  empirical 
description  of  the  data  set.  Under  no  lack-of-fit  conditions,  the  mean  square 
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error  and  mean  square  pure  error  should  doth  estimate  the  population  variance 
a2.  Inaeed,  MS p  =  4452.5  and  MS  =  4443.1*  producing  estimated  sample  standard 
deviations  of  56.30  pg/mL  and  66.66  pg/mL.  These  standard  deviations  were 
somewhat  larger  than  those  listed  in  tne  literature  (by  the  criterion  that  a 
normal  range  covers  a  oistance  of  accut  4o,  the  stanaard  deviation  incicatea  by 
the  literature  is  about  20  pg/mL).  -’owever,  the  given  estimates  of  :  were 
inflated  by  the  presence  of  potential  outliers.  Since  the  value  of  Cock's  C  was 
not  considered  extreme  (all  had  Cook's  Os  of  between  0.01  and  0.04)*  the  4 
possible  outliers  (corresponding  to  animal  130  (week  -3).  animal  153  (week  -3), 
animal  159  (week  -2)*  and  animal  134  (weex  O)  were  not  rejectee  from  the  data 
set.  The  high  values  of  these  ccservatlons  (all  above  600  pg/mL)  did  tend  to 
raise  the  mean  values  at  those  weeks*  and  thus  inflated  the  estimates  of  the 
standard  deviation. 

The  next  model  verification  step  involved  examining  the  residual  and 
partial  residual  plots  to  confirm  the  least  squares  regression  assumption  that 
the  model  errors  were  NID(0*<?2).  This  step  wculc  defend  the  use  of  F  tests  to 
determine  tne  statistical  significance  of  the  parameters.  Additionally*  this 
step  would  validate  the  statistics  which  produced  tables  listing  confidence 
intervals  of  the  norepinephrine  concentrations.  A  number  of  residual  plots 
suggested  themselves  immediately:  residuals  versus  time,  residuals  versus 

predicted  value  of  norepinephrine  concentration,  residuals  versus  animal  case 
number*  studentized  residuals  versus  the  previous  three*  and  partial  residual 
plots  corrected  for  the  model  terms  3g,  3^,  and  3  Examination  of  the 

residual  plots  yielded  no  discernible  patterns  in  the  distribution  of  the 
residuals.  Thus,  the  resiauals  were  normally  distributed  with  mean  0  and 
variance  a2.  The  residual  plots  are  included  in  Appendix  F. 

Since  the  data  from  this  study  arose  as  a  time  series*  there  was  a 
possibility  that  the  residuals  were  in  some  part  autocorrel ated  to  prior 
observations.  To  determine  the  extent  of  this  autocorrelation,  an 
autoregressive  model  building  procedure  (PR0C  AUTOREG  from  the  SAS  ETS  series) 
was  used  with  lag  times  of  0,  1,  2*  3,  and  4  weeks. 

Results  of  the  autoregression  (not  included  in  report)  indicated  a 
significant  amount  of  correlation  between  data  at  one  week  to  data  at  the 
previous  week  (lag-1  autocorrelation).  The  autoregression  also  detected  a 
considerably  smaller  (although  statistically  significant)  lag-2  autocorrelation. 
Lag-1  correlation  indicated  that  the  best  predictor  of  any  single  observation 
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was  the  previous  observation  for  that  particular  animal  (rather  than  the  value 
yielced  by  substituting  the  parameter  estimates  and  week  number  into  the  derived 
norepinephrine  model).  If  the  study  purpose  were  to  determine  a  predictive 
model  cf  norepinephrine  behavior  in  the  rats,  then  the  just  mentioned  conclusion 
would  have  dire  consequences  with  regards  to  the  model  obtained  in  Appendix  D. 
However,  the  main  reason  regression  was  chosen  to  model  this  data  was  not  to 
produce  a  precictive  model  of  norepinephrine  versus  time,  but  rather  to 
determine  whether  or  not  two  blocked  groups  (exposure  and  sham-exposure) 
cisplayed  any  differences  in  norepinephrine  behavior.  For  this  purpose,  non¬ 
independence  in  the  residuals  does  not  call  into  question  the  overall 
conclusions  drawn  from  the  model.  To  compensate  for  this  deficiency,  it  would 
only  be  necessary  to  raise  the  ct  risk  used  in  determining  the  norepinephrine 
model.  Since  3g,  B 1  ,  and  3^  were  found  to  be  significant  at  probabilities 
less  than  0.0001,  this  alteration  of  significance  had  no  practical  effect  on  the 
final  model  determined  in  the  analysis.  The  large  number  of  observations  taken 
essentially  made  this  data  set  relatively  insensitive  to  potential  problems 
(such  as  lack-of-fit  or  non  independence  in  the  residuals). 

To  complete  the  analysis,  diagnostics  to  check  for  model  multicoll inearity 
and  correlation  between  the  terms  were  used.  Examination  of  the  listed 
condition  numbers  and  matrix  eigenvalues  (being  provided  under  separate  cover) 
detected  no  troublesome  values.  This  review  indicated  that  the  model  did  not 
display  a  significant  degree  of  multicoll inearity .  Similarly,  examination  of 
the  correlation  matrix  showed  that  correl ations between  the  estimated  values  cfS 
were  all  within  tolerable  limits.  The  highest  degree  of  correlation  was  between 
the  x  and  the  x  term,  which  often  occurs  when  using  a  polynomial  model  In 
linear  regression. 

^or  future  reference,  and  for  the  sake  of  completeness,  tables  listing 
animal  case  number,  observations  (if  taken)  at  each  week,  predicted  value  of 
norepinephrine  concentration,  standardized  error  of  prediction,  95 %  confidence 
intervals  on  the  mean  value  of  the  norepinephrine  concentration,  and  residuals 
were  prepared,  as  were  tables  containing  animal  case  number,  regular  and 
studentized  residual  values,  a  graphical  display  of  student  residual  values,  and 
influence  statistics  (such  as  Cook's  0).  These  tables  were  used  to  detect  both 
outliers  and  influential  data  points  in  the  norepinephrine  data  set. 

To  arrive  at  a  conservative  estimate  of  the  minimum  change  due  to  RFR  in 
resting  norep inepnrine  concentrations  which  this  protocol  was  capable  of 


detecting,  the  value  of  the  operating  curve  parameter  Jg  corresponding  to  the 
RFR  factor  (B)  discussed  at  the  beginning  of  the  statistical  methodology  was 
calculated.  This  parameter  was  given  by 


,2  naD2 

Q  S  ■  —  — — ■ 

B  2bc2 


(B-7) 


where  n  =  number  of  replications  per  cell  =  40, 
a  =  number  of  levels  of  factor  A  =  12, 

b  =  number  of  levels  of  factor  B  =  2, 

a 2  *  population  variance,  and 

=  detection  threshold. 

Substituting  In  values  fcr  a,  b,  n,  and  the  MS 
provided  an  operating  curve  parameter  of 


pe 


n 

as  an  estimate  of  a 


*B  =  0.1643  D.  (B-8) 

To  obtain  a  value  of  5  from  the  operating  curve,  the  type  I  risk  a  and  type  II 
risk  3  were  set  to  0.05  and  0.10,  respectively.  Then,  the  value  $  was  read  from 
the  fixed  effects  ANOVA  curve  with  ■  1  and  Vj  *  936.  This  value  was 

3  2.4.  (B-9) 

Note  that  the  degrees  of  freedom  for  the  numerator,  and  the  degrees  of 
freedom  for  the  denominator,  Vj,  #ere  calculated  with  the  equation 

^  »  b-1,  and 
^2  =  ab(n-l). 

The  detection  level  was  therefore 

Dg  =  14.60  pg/ml.  (B-12) 

Thus,  this  protocol  was  able  to  conservatively  detect  an  Increase  In 
resting  plasma  norepinephrine  concentrations  of  14.60  pg/ml  about  905!  of  the 
time. 


(B-10) 

(B-ll) 
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In  many  ways/  the  epinephrine  scatter  diagrams  of  Figure  7  ar.d  8  closely 
resemoled  the  norepinephrine  scatter  diagrams.  Therefore#  epinephrine 
concentration  was  modelled  In  a  similar  manner  to  norepinephrine.  The  equation 


y  =  +  «qZ  +  3^zx  +  a^zx2  C B— 13 ) 

where  y  =  plasma  epinephrine  concentration  (In  pg/mL) 
x  =  time  (In  weeks)#  and 

z  =  a  categorical  variable  with  value  0  for  animals  In  the  sham- 
exposure  group  and  value  I  for  animals  In  the  exposure  group# 

was  tested  for  the  significance  of  the  coefficients  ag»  a^»  and  a^.  These 
terms  described  the  RFR- Interact Ion  with  the  resting  epinephrine  concentration. 

Data  from  the  epinephrine  spreadsheets  (Appendix  G)  were  subsequently  put 
Into  a  new  file  and  a  second  SAS  formatting  program  (included  In  Appendix  H)  was 
prepared  to  analyze  the  data. 

The  model  Indicated  by  the  forward  stepwise  and  maximum  R^  regression 
procedures  was 


where 


B0  +  ®lx 


3  x2, 
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S 

II 


158.80, 
-6.62#  and 
0.28, 


(B-14) 


with  the  x#  y#  and  z  variables  as  defined  previously.  The  entry  and  exit  risk 
were  both  set  to  0.095.  The  outputs  of  both  regression  procedures  are  Included 
In  Appendix  I.  Note  that  the  absence  of  a  terms  Indicated  that#  at  a  risk  of 
0.095#  there  was  no  statistical  difference  In  plasma  epinephrine  concentrations 
between  the  exposure  and  sham-exposure  groups.  The  estimated  resting 
concentration  of  plasma  epinephrine,  158.8  pg/mL,  also  agreed  well  with 
established  values  cited  In  the  literature  (180  +35  pg/mL).  This  agreement  was 
a  further  Indication  of  no  systematic  error  within  the  sampling/assaying 
procedure. 

Epinephrine  concentration  In  the  sham-exposure  and  exposure  groups 
displayed  the  same  type  of  time  dependency  found  In  the  norepinephrine 
concentrations.  Since  epinephrine  and  norepinephrine  release  within  the  DOdy 
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are  physiologically  coupled,  this  was  rot  a  surprising  find.  Specifically, 
resting  epinephrine  values  were  at  their  highest  value  of  181.3  pg/mL  at  the 
study  onset  (week  -3).  The  resting  level  then  gradually  declined,  reaching  its 
lowest  point  of  119.4  pg/ml  at  week  12  of  the  study.  Epinephrine  concentration 
slowly  rose  beyond  that  point  to  a  value  of  133.9  pg/ml  at  week  19,  the  last 
week  for  which  data  were  taken.  All  of  the  just  mentioned  values  are  typical  of 
resting  epinephrine  concentrations  in  normal,  unstressed  rats.  It  did  not 
appear,  therefore,  that  chronic  exposure  to  435-MHz  RFR  Induced  any  stress,  as 
measured  by  the  resting  concentration  of  plasma  epinephrine,  in  the  exposure 
group  when  compared  to  the  sham-exposure  group. 

The  just  mentioned  conclusions  could  only  be  accepted  upon  verification  of 
the  assumptions  used  in  building  the  model.  These  assumptions  Included  no  model 
lack-of-fit,  NID(0,CT^)  residual  distribution,  and  no  model  multicoll Inearity. 

First,  the  model  was  checked  for  lack-of-fit  (Appendix  J).  The  mean  square 
error  and  mean  square  pure  error  were  3359.31  and  3296.69  respectively,  yielding 
sample  standard  deviation  estimates  of  57.96  pgVmL  and  57.42  pg/mL.  Since  ooth 
of  these  estimates  were  rather  close  to  one  another,  lack-of-fit  was  probably 
not  significant.  The  computed  lack-of-fit  test  statistic  was  then  found  to  be 
smaller  than  the  critical  value.  This  test  confirmed  that  model  lack-of-fit  was 
not  present. 

The  epinephrine  data  set  was  then  checked  for  outlier  data  values  before 

generating  residual  plots.  Three  observations  at  week  -3  (animal  #53, 

[epinephrine]  =  560  pg/mL;  animal  #57,  [epinephrine]  =  806  pg/mL;  and  animal 

#62,  [epinephrine]  *  540  pg/mL)  were  determined  to  be  outliers  and  were 

subsequently  removed  from  the  data  set.  All  three  points  had  values  of  Cook's  D 

greater  than  0.05,  and  thus  were  overly  Influential  In  comparison  with  other 

data  points  from  week  -3.  Once  the  data  set  was  edited,  residual  plots  were 

generated  to  check  the  assumption  that  the  model  errors  were  distributed  NID(0» 

o  ) .  Appendix  K  contains  the  epinephrine  residual  plots.  Examination  of  the 

plots  yielded  no  obvious  patterns  or  problems,  thereby  Indicating  that  the 

2 

residuals  were  normally  distributed  with  mean  0  and  variance  a  . 

However,  independence  among  the  residuals  was  not  assured.  Often, 
residuals  produced  from  a  regression  modelling  data  taken  in  time  series  show  a 
degree  of  autocorrelation  from  one  week  to  the  next.  To  adequately  address  this 
problem,  it  became  necessary  to  perform  an  autoregression  on  the  regression 


model,  and  determine  the  extent  of  autocorrelation  and  the  effects  of  the 
autocorrelation  on  the  hypothesis  tests. 

Results  of  the  autoreg ress 1  on  (not  included  in  text)  indicated  a 
significant  amount  of  correlation  between  data  at  one  week  to  data  at  the 
previous  week  (lag-1  autocorrelation).  The  autoregression  also  detected  a 
smaller  amount  of  lag-3  correlation  stemming  from  a  presently  unknown  source. 
The  1  ag-1 '  correl  ation  Indicated  that  the  best  predictor  for  an  animal’s 
epinephrine  concentration  was  more  the  last  known  epinephrine  concentration 
rather  than  the  time  into  the  study.  If  the  study  purpose  were  to  determine  a 
predictive  model  of  epinephrine  concentration  versus  time,  this  would  be  a 
significant  find.  However,  since  the  study  purpose  was  to  determine  the  effects 
of  RFR  on  epinephrine  concentration,  this  finding  did  not  significantly  change 
any  of  the  significance  tests  on  the  parameters.  To  adjust  for  the  presence  cf 
lag-1-  correlation  on  the  parameter  tests,  a  qualitative  measure  would  be  to 
increase  the  *  risk  of  the  conclusion  drawn.  Since  the  probability  that  each 
epinephrine  parameter's  F  statistic  is  greater  than  Fq  was  better  than  0.00C1, 
then  this  adjustment  of  a  risk  would  have  no  practical  effect  on  the  conclusion. 
Thus,  although  the  model  had  nonindependent  characteristics,  they  were  of  such  a 
nature  as  to  not  affect  the  final  conclusion  taken  from  the  model. 

To  complete  the  analysis,  diagnostics  to  check  for  model  multicoll inearity 

and  correlation  between  the  terms  were  used.  Examination  of  the  listed 

condition  numbers  and  matrix  eigenvalues  (being  provided  under  separate  cover) 

detected  no  troublesome  values.  This  review  Indicated  that  the  model  did  not 

display  a  significant  degree  of  multicoll Inearity.  Similarly,  examination  of 

the  correlation  matrix  showed  that  correlation  between  the  estimated  values  of  3 

were  all  within  tolerable  limits.  The  highest  degree  of  correlation  was  between 
2 

the  x  and  the  x  term,  which  often  occurs  when  using  a  polynomial  model  in 
linear  regression. 

For  future  reference,  and  for  the  sake  of  completeness,  tables  listing 
animal  case  number,  observations  (If  taken)  at  each  week,  predicted  value  of 
epinephrine  concentration,  standardized  error  of  prediction,  95 %  confidence 
intervals  on  the  mean  value  of  the  epinephrine  concentration,  and  residuals  were 
prepared,  as  were  tables  containing  animal  case  number,  regular  and  studentized 
residual  values,  a  graphical  display  of  student  residual  values,  and  influence 
statistics  (such  as  Cook's  D).  These  tables  were  used  to  detect  both  outliers 
and  influential  data  points  in  the  epinephrine  data  set. 
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To  arrive  at  a  conservative  estimate  of  the  minimum  cnance  due  to  RFR  in 
resting  epinephrine  concentrations  which  this  protocol  was  capable  of  detecting* 
the  value  of  the  operating  curve  parameter  corresponding  to  the  RFR  factor 
(3)  discussed  at  the  beginning  of  the  statistical  methodology  was  calculated. 
This  parameter  was  given  by 


naD 

2bo" 


where  n  =  number  of  replications  per  cell  =  40, 
a  =  number  of  levels  of  factor  A  =  12, 


(B-1S) 


b  =  number  of  levels  of  factor  8=2* 
a  =  population  variance,  and 
=  detection  thresnola. 


Suostituxing  in  values  for  a,  o,  n, 
proviGec  an  operating  curve  parameter  of 


and  the  MS  as  an 
pe 


estimate  of 


cr2  . 


*B  *  0.1908  D.  (3-16) 

To  obtain  a  value  of  ?  from  the  operating  curve,  the  type  I  risk  a  and  type  II 
risk  3  were  set  to  0.C5  and  0.10,  respectively.  Then,  the  value  of  $  was  read 
from  the  fixed  effects  ANOVA  curve  with  =  l  and  ^  =  936.  This  value  was 


*s  =  2.4.  f  B— 17 ) 

Degrees  of  freedom  in  both  numerator  and  denominator  were  calculated  in  the 
same  manner  as  those  in  the  norepinephrine  analysis.  Note  that  the  40 
replications  in  the  protocol  were  not  replications  in  the  truest  sense  of  the 
word  (since  a  single  animal  was  not  put  through  the  study  40  times).  Since 
Sprague-Dawley  rats  represented  a  very  homogeneous  population,  this  difference 
would  have  only  minor  effects  on  the  rigor  of  this  calculation. 

The  detection  level  was  therefore 


Dg  =  12.58  pg/mt. 


(B-1S) 
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Thus/  this  protocol  conservatively  was  able  to  detect  an  increase  in 
resting  plasma  epinephrine  concentrations  of  12.58  pg/mL  about  90%  of  the  time. 

Plasma  Oooamine  Statistical  Analysis. 

Upon  examining  the  scatter  diagrams  of  Figures  10  and  11  and  the  mean 
dopamine  concentration  versus  time  plot  of  Figure  12/  it  did  not  appear  that 
resting  dopamine  levels  in  the  exposure  group  were  higher  than  resting  dopamine 
levels  in  the  sham-exposure  group.  Therefore/  the  model  to  test  for  RFR-induced 
effects  on  dopamine  concentration  was  the  starting  model  of  the  norepinephrine 
and  epinephrine  analyses: 


y 


+ 


(5-19) 


where  y  =  resting  plasma  dopamine  concentration  (in  pg/rnL), 
x  =  time  (in  weeks)/  and 

z  =  a  categorical  variable  with  value  0  for  animals  in  the  sham- 
exposure  group  ana  value  1  for  animals  in  the  exposure  group. 


The  significance  of  the  a  terms  in  this  model  determined  whether  or  not  there 
were  any  RFR-induced  effects;  the  algebraic  sign  of  the  a  then  oexermined 
whether  or  not  the  effects  tended  to  increase  resting  hormone  concentrations 
(Indicated  by  positive  a)  or  decrease  resting  hormone  concentrations  (Indicated 
by  negative  a).  Note  that  these  ot  terms  should  not  be  confused  with  the  symbol 
for  statistical  significance  (risk)/  which  is  also  ana. 

Data  from  the  dopamine  spreadsheets  (Appendix  L)  Were  subsequently  put  into 
a  new  file  and  a  tnira  SAS  formatting  program  (Appendix  M)  was  prepared  to 
analyze  the  data. 

The  model  Indicated  by  the  forward  stepwise  and  maximum  R  regression 
proceuures  was 
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where 


S0  =  51.19, 

3*  =  -3.14, 

X 

*i  *  -0.92,  anc 
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with  the  x,  y,  and  z  variables  defined  as  previously.  The  entry  and  exit  risk 
were  both  set  to  0.10.  The  outputs  of  both  regression  procedures  are  included 
in  Appendix  N.  The  absence  of  indicated  that  RFR  did  not  produce  a 
detectable  effect  on  the  intercept  of  the  model*  ana  therefore  did  not  bias  the 
dopamine  concentration  of  the  exposure  group  when  compared  to  the  snam-exposure 
group.  Equivalently*  this  showed  that  at  the  onset  of  exposure  (week  0)*  both 
groups  displayed  comparable  resting  dopamine  levels.  This  result  was  not 
surprising*  since  tne  experiment  was  designed  such  that  the  initial  resting 
copamine  levels  of  both  groups  wculc  be  similar.  Additionally,  there  w as  no 
evidence  of  any  RFR-induced  effect  on  the  curvature  of  the  exposure  group. 

The  exposure  group  did  differ  from  the  sham-exposure  group  with  regards  to 
overall  time  response,  however.  In  both  groups*  dopamine  concentration  started 
out  somewhat  high  (61.8  pg/mL  sham-exposure  group;  64.6  pg/mL  exposure  group  at 
week  -3).  After  the  initiation  of  radiation,  the  exposure  group's  estimated 
resting  dccamine  concentration  remained  below  that  of  the  sham-exposure  group 
for  the  duration  of  the  study.  At  week  12*  estimated  resting  dopamine 
concentrations  in  sham-exposure  animals  reached  a  low  value  of  32.3  pg/mL;  tne 
low  for  exposure  animals  was  attained  at  week  16  with  an  estimated  dopamine 
level  of  15.6  pg/mL.  The  dopamine  concentration  then  rose  slightly,  reaching 
estimated  values  of  35.6  pg/mL  in  snam-exposure  animals  ana  21.1  pg/'niL  in 
exposure  animals  oy  week  19  (the  final  week  data  were  collected)  of  the  study. 
Both  ranges  (32.3  to  61.8  pg/mt  in  sham-exposure  animals*  19.6  to  64.6  pg/mL  in 
exposure  animals)  were  still  well  within  the  normal  range  of  plasma  dopamine  in 
ncnstressed  male  Sprague-Dawley  rats  (85  +35  pg/mL).  Stress  in  these  animals 
is  reflected  in  an  Increased  rate  of  dopamine  secretion.  Therefore,  these 
results  indicated  that  chronic  exposure  to  435-MHz  RFR  did  not  induce  an 
elevati-cn  in  resting  dopamine  concentration  in  the  exposure  group. 

Once  again,  it  was  then  necessary  to  check  the  validity  of  the  assumptions 
used  in  building  the  dopamine  regression.  First*  a  model  lack-of-fit  test  was 
performed  (Appendix  0).  The  mean  square  error  and  the  mean  square  pure  error 
were  1235.37  and  314.01  respectively,  yielding  sample  stanaard  deviation 
estimates  of  35.15  and  28.53  pg/mL.  The  calculated  value  of  Fg  was  then  about 
1.51,  wnile  the  critical  value  was  about  1.38. 


Since  Fg  exceeded  the  critical  value,  the  dopamine  model  displayed  a 
significant  lack-of-fit,  thereby  deviating  from  results  obtained  in  plasma 
noreo’ nephrine  and  plasma  epinephrine.  The  situation  was  reminiscent  of  that 
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encountered  in  the  analysis  of  ACTH  and  corticosterone/  and  in  the  analysis  of 
prolactin  C2»3 ].  In  those  cases/  significant  lack-of-fit  was  handled  by 
qualitatively  altering  the  significance  levels  ct  to  compensate  for  the  model 
defects.  This  procedure  was  preferable  to  transformation  of  the  dependent  or 
independent  variables/  since  a  transformation  on  the  dependent  variable  y  would 
alter  the  residual  distribution  and  a  transformation  on  the  independent  variable 
x/  although  theoretically  possible/  would  be  time  consuming  and  costly  and  yield 
a  model  with  minimally  better  predictive  value. 

We  then  decided  to  follow  this  course  for  the  dopamine  model.  Therefore/ 
model  lack-of-fit  could  be  deemed  statistically  significant  but  practically 
insignificant  by  altering  the  a  risk  in  the  coefficients.  Since  those  co¬ 
efficients  were  highly  significant  to  begin  with/  this  alteration  of  a  risk 
should  not  change  the  model  in  any  manner. 

Residual  plots  were  then  generated  for  the  dopamine  data.  Since  no 
observations  in  the  original  data  set  had  values  of  Cook’s  D  higher  than  0.05/ 
we  decided  not  to  reject  any  values  from  the  data  set.  The  residual  plots 
(Appendix  P)  therefore  displayed  no  obvious  patterns  or  problems.  This 
supported  the  assumption  that  the  model  errors  were  normally  distributed  with  a 
mean  of  zero  and  a  variance  of  cr  . 

As  previously  mentioned/  the  lack  of  patterns  within  the  residual  plots  did 
not  guarantee  Independence  within  the  observations  because  models  produced  by 
regression  of  data  taken  in  time  series  tend  to  show  some  degree  of 
autocorrelation  between  the  £^s  of  each  time  Interval.  To  adequately  address 
this  question/  the  dopamine  data  set  was  reexamined  with  an  autoregressive 
procedure  to  determine  the  extent  of  residual  autocorrelation  and  its  effects  on 
the  model's  hypothesis  tests. 

Results  of  the  autoregression  (not  included  in  this  report)  Indicated  a 
significant  amount  of  correlation  within  the  data  at  lags  of  1  and  2  weeks  (the 
week  2  autocorrelation  was  considerably  smaller  than  the  week  1  autocorrelation/ 
and  stems  from  a  presently  unknown  source).  Once  again/  this  quantitative 
estimate  of  autocorrelation  was  not  unexpected/  nor  practically  significant  in 
terms  of  the  conclusions  drawn  from  the  model.  A  further  adjustment  of  the  a 
risk  values  in  the  regression  would  compensate  for  the  lag-1  autocorrelation. 
Since  the  probability  that  each  parameter  F  statistic  was  greater  than  the 
critical  F  value  was  better  than  0.0003  (for  the  parameters  statistically 
significant  in  the  dopamine  regression)/  this  adjustment  of  -  risk  was 
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inconsequential.  Thus*  the  sheer  number  of  observations  taken  helped  compensate 
for  the  model's  two  main  defects:  lack-of-fit  and  nonindependent  residuals. 

To  complete  the  analysis*  diagnostics  to  check  for  model  multicoll inearity 
and  correlation  between  the  terms  were  used.  Examination  of  the  listed 
condition  numbers  and  matrix  eigenvalues  (being  provided  under  separate  cover) 
detected  no  troublesome  values  and  Indicated  that  the  model  did  not  display  a 
significant  degree  of  mul t Icoll 1  near  1 ty .  Similarly*  examination  of  the 
correlation  matrix  showed  that  correl ations  between  the  estimated  values  of  6 
were  all  within  tolerable  limits.  The  highest  degree  of  correlation  was  between 
the  x  and  the  x^  term*  which  often  occurs  when  using  a  polynomial  model  in 
linear  regression. 

For  future  reference*  and  for  the  sake  of  completeness*  tables  listing 
animal  case  number,  observations  (if  taken)  at  each  week,  predicted  value  of 
dopamine  concentratl on*  standardized  error  of  prediction,  95%  confidence 
intervals  on  the  mean  value  of  the  dopamine  concentration,  and  residuals  were 
prepared*  as  were  tables  containing  animal  case  number,  regular  and  studentized 
residual  values,  a  graphical  display  of  student  residual  values,  and  Influence 
statistics  (such  as  Cook's  D),  These  tables  were  used  to  detect  both  outliers 
and  influential  data  points  in  the  dopamine  data  set. 

To  arrive  at  a  conservative  estimate  of  the  minimum  change  due  to  RFR  in 
resting  dopamine  concentrations  which  this  protocol  was  capable  of  detecting* 
the  value  of  the  operating  curve  parameter  corresponding  to  the  RFR  factor 
(B)  discussed  at  the  beginning  of  the  statistical  methodology  was  calculated. 
This  parameter  was  given  by 

.  _2 
.  2  naD 

’b  '  "ST"  !B-21) 

where  n  -  number  of  replications  per  cell  =  40, 
a  =  number  of  levels  of  factor  A  =  12* 
b  =  number  of  levels  of  factor  B  *  2* 

c  =  population  variance,  and 

2 

0*  =  detection  threshold. 

Substituting  in  values  for  a,  b*  n,  and  the  MSpe  as  an  estimate  of  ^ 
provided  an  operating  curve  parameter  of 
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ig  =  0.3840  D. 


(3-22) 


To  obtain  a  value  of  f  from  the  operating  curve*  the  type  I  risk  a  and  type  II 
risk  3  were  set  to  0.05  and  0.10,  respectively.  Then,  the  value  of  3  was  read 
from  the  fixed  effects  ANOVA  curve  with  =  l  and  ^  =  936.  This  value  was 

5g  =  2.4.  (3-23) 

Degrees  of  freedom  in  both  numerator  and  denominator  were  calculated  in  the 
same  manner  as  those  in  the  norepinephrine  analysis.  Once  again,  the  40 
replications  in  the  protocol  were  not  renl ications  in  the  truest  sense  (since  an 
individual  animal  was  not  put  through  the  study  40  times).  However,  Sprague- 
Dawley  rats  represent  a  very  homogeneous  population  and  thus  minimize  the 
between-individual  variation  of  the  cell  observations. 

The  detection  level  was  therefore 

Dg  =  6.25  pg/ml.  (B-24) 

Thus,  the  protocol  was  able  to  detect  an  increase  in  resting  plasma 
dopamine  concentrations  of  6,25  pg/ml  about  90%  of  the  time. 

We  gratefully  acknowledge  the  assistance  of  Dr.  Russell  G.  Heikes  of 
Georgia  Tech's  Department  of  Industrial  and  Systems  Engineering  in  developing 
the  statistical  methodology  of  this  appendix. 
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NAME  ME.MTYPE  OBS  TRACKS  PROT 

OVLMN  DATA  43  i 

SCTR  Data  3024  i 
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3  SAS(R)  LOG  CMS  SAS  5.16  VM/CMS  CMS  USER  OSECLSB 


SNORE? IN/ DATA 

3072 

1 

TESTN  /DATA 

3024 

1 

72  DELETE 

’3  DELETE 

LIST  OF  MEMBERS 

SCTR ; 

OVLMN; 

AFTER  UPDATE  OF 

DIRECTORY. 

NAME  MEMTYPE 

OBS 

TRACKS  PROT 

SNOREPIN/DATA 

3072 

1 

TESTN  /DATA 

3024 

1 

NOTE:  THE  PROCEDURE  DATASETS  USED  0.12  SECONDS  AND  424K. 

PROC  STEPWISE: 

75  MODEL  Y  -  X  XSQR  2  XZ  XSQRZ  /SLENTR Y-0 . 10  SLSTAY=0.10  STEPWISE  MAXR; 
NOTE:  THE  PROCEDURE  STEPWISE  USED  0.63  SECONDS  AND  424K  AND  PRINTED  PAGES  6  TO  8. 

76  PROC  PRINTTO  NEW  UNIT-23; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  .AND  424K. 

77  PROC  REG; 

78  MODEL  Y  -  X  XSQR  /  PARTIAL; 

79  ID  CASE; 

NOTE:  ACOV  AND  SPEC  OPTION  ONLY  VALID  WITH  RAWDATA 

NOTE:  THE  PROCEDURE  REG  USED  1.46  SECONDS  AND  744K  AND  PRINTED  PAGES  9  TO  12. 

80  PROC  PRINTTO  NEW  UNIT-24; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  424K. 

81  PROC  GLM; 

82  CLASS  X  Z; 

83  MODEL  Y  -  X  X*X  X*Z; 

NOTE:  THE  PROCEDURE  GLM  USED  3.36  SECONDS  AND  1128K  AND  PRINTED  PAGES  13  TO  14. 

84  PROC  PRINTTO  NEW  UNIT-25; 

85  K - * 

86  *  * 

87  -  to  obtain  tables  listing  the  variance  inflation  factors,  * 

88  *  influence  statistics,  and  tolerances,  the  following  SAS  * 

89  *  statements  were  used  in  this  partition:  * 

90  * 

91  *  PROC  REG; 


92  *  MODEL  Y  -  X  XSQR  /  TOL  VIF  INFLUENCE;  * 

>3  *  ID  CASE;  * 

*  OUTPUT  OUT-RNOREPIN  P-PREDICT  R-RESID  STUDENT-STUDENT ; 

95  * 

95  Dr - *. 


NOTE:  THE  PROCEDURE  PRINTTO  USED  0.04  SECONDS  AND  424K. 

97  PROC  REG; 

98  MODEL  Y  -  X  XSQR  /  I  SSI  SS2  STB  COVB  CORRB  SEQB  COLLIN 

99  COLLINOINT  ACOV  P  R  CLM; 

100  ID  CASE; 

101  OUTPUT  OUT-RNOREPIN  P-PREDICT  R-RESID  STUDENT-STUDENT; 

NOTE:  THE  DATA  SET  WORK. RNOREPIN  HAS  3072  OBSERVATIONS  AND  13  VARIABLES. 

NOTE:  THE  PROCEDURE  REG  USED  6.76  SECONDS  AND  744K  .AND  PRINTED  PAGES  15  TO  80. 

132  PROC  PRINTTO  NEW  UNIT-26: 
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SAS(R)  LOG  CMS  SAS  5.16 
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NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  424K. 

103  PROC  PLOT  DATA-RNOREPIN; 

104  PLOT  RESID*X-’*’  /  V AXIS— 200  TO  200  BY  50; 

105  PLOT  RESID”PREDICT= ' ” '  /  HAXIS-22G  TO  300  BY  5  VAXIS— 200  TO  200  BY  50; 

106  PLOT  STUDENT *X« ’  /  VAXIS— 3  TO  3  BY  0.5; 

107  PLOT  STUDENT “PREDICT- ' * '  /  HAXIS-220  TO  300  BY  5  VAXIS— 3  TO  3  BY  0.5; 

108  TITLE  'NOREPINEPHRINE  RESIDUAL  PLOTS’; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.34  SECONDS  AND  424K  AND  PRINTED  PAGES  81  TO  84. 

109  PROC  PRINTTO  NEW  UNIT-27 ; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  424K. 

110  PROC  PLOT  DATA-RNOREPIN; 

111  BY  2; 

112  PLOT  RESID*CASE-'*'  /  VAXIS— 200  TO  200  BY  50  HAXIS-0  TO  65  BY  5; 

113  PLOT  STUDENT*CASE-’*'  /  VAXIS— 3  TO  3  BY  0.5  HAXIS-0  TO  65  BY  5; 

114  TITLE  ’NOREPINEPHRINE  RESIDUAL  PLOTS’; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.74  SECONDS  AND  424K  AND  PRINTED  PAGES  85  TO  88. 

115  PROC  PRINTTO  NEW  UNIT-28 ; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  424K. 

116  PROC  AUTOREG; 

117  TITLE  'Norepinephrine  Autoregressive  Models'; 

118  MODEL  Y  -  X  XSQR  /  COEF  CORRB  COVB  BACKSTEP; 

119  MODEL  Y  -  X  XSQR  /  NLAG-1  COEF  CORRB  COVB  BACKSTEP; 

120  MODEL  Y  -  X  XSQR  /  NLAG-2  COEF  CORRB  COVB  BACKSTEP; 

121  MODEL  Y  -  X  XSQR  /  NLAG-3  COEF  CORRB  COVB  BACKSTEP; 

122  MODEL  Y  -  X  XSQR  /  NLAG-4  COEF  CORRB  COVB  BACKSTEP; 

NOTE:  THE  PROCEDURE  AUTOREG  USED  6.64  SECONDS  AND  424K  AND  PRINTED  PAGES  89  TO  101. 
NOTE:  SAS  USED  6952K  MEMORY. 

NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 

PO  BOX  8000 

CARY,  N.C.  27511-8000 
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THE  ABOVE  MODEL  IS  THE  BEST  5  VARIABLE  MODEL  FOUND. 


APPENDIX  E 

NOREPINEPHRINE  LACK-OF-FIT  TEST 


APPENDIX  F 

NOREPINEPHRINE  RESIDUAL  PLOTS 
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:  SaS(R)  LOG  CMS  SAS  5  - 16  VM/ CMS  CMS  USER  QSECLSB 

NOTE:  COPYRIGHT  > C)  1984,1986  SAS  INSTITUTE  INC.,  CARY,  N.C.  27511,  U.S.A. 
NOTE:  CMS  SaS  RELEASE  5.16  AT  GEORGIA  INSTITUTE  OF  TECHNOLOGY  (03559001). 

NOTE:  CPUID  VERSION  »  FF  SERIAL  -  012242  MODEL  -  4381  . 

NOTE:  SAS  OPTIONS  SPECIFIED  .ARE: 

LEAVE-0 

1  data  teste; 

2  CMS  FILEDEF  X  DISK  EPIN  DAT  Ai ; 


3 

CMS 

FILEDEF 

20 

DISK 

EPINO 

LISTING 

Al 

( 

3LKSIZE 

141 

RECFM 

VB.A 

LRECL 

133 

4 

CMS 

FILEDEF 

21 

DISK 

EPIN1 

LISTING 

Al 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

5 

CMS 

FILEDEF 

22 

DISK 

EPIN2 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

6 

CMS 

FILEDEF 

23 

DISK 

EPIN3 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

CMS 

FILEDEF 

24 

DISK 

EPIN4 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

8 

CMS 

FILEDEF 

25 

DISK 

EPIN5 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

9 

CMS 

FILEDEF 

26 

DISK 

EPIN6 

LISTING 

Al 

( 

\ 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

10 

CMS 

FILEDEF 

27 

DISK 

EPIN  7 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

11 

CMS 

FILEDEF 

28 

DISK 

EPIN8 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

12  .ARRAY  WEEK  (24}  WKN3  WKN2  MISSN1  WK0-WK20; 

13  KEEP  X  XSQR  Y  2  XZ  XSQR2  CASE: 


1- 

INFILE  X: 

INPUT  CASE 

1-3 

16 

WKN3 

5-7 

i  / 

WKN2 

9-11 

18 

WKO 

13—15 

19 

WK1 

17-19 

20 

WK2 

21-23 

21 

WK3 

25-27 

22 

WK4 

29-31 

23 

WK5 

33-35 

24 

WK6 

37-39 

25 

WK7 

41-43 

26 

WK8 

45-47 

n  - 

/ 

WK9 

49-51 

28 

WK10 

53-55 

29 

WK1 1 

57-59 

30 

WK12 

61-63 

31 

WK13 

65-67 

32 

WK14 

69-71 

33 

WK15 

73-75 

34 

WK16 

77-79 

35 

WK1 7 

81-83 

36 

WK18 

85-87 

3' 

WK19 

89-91 

38 

WK20 

93-95 

39  ; 

**0  MISSN1-.  ; 

-1  IF  CASE  <  100  THEN  2-0; 

42  IF  CASE  >-  IOC  THEN  Z  -  1; 

“3  IF  2-1  THEN  CaSE-CaSE- 100 ; 

44  DC  I  -  1  TO  24; 

-*5  X  -  1-4;  XSQR  -  X*X;  XZ  -  X*Z;  XSQRZ  -  X,rX*Z;  Y  -  WEEK  111  ; OUTPUT; 
46  END; 


S  r:-E  E ?  . N  0 A 7 

WERE  RE.-uD  FROM  INFILE  X 


30 


SAS(R)  LOG  CMS  SAS  5.16 


VM/ CMS  CMS  USER  QSECLSB 


NOTE:  DATA  SET  WORK. TESTE  HAS  3096  OBSERVATIONS  .AND  7  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.59  SECONDS  -AND  208K. 

47  ?ROC  CONTENTS; 

NOTE:  THE  PROCEDURE  CONTENTS  USED  0.20  SECONDS  .AND  464K  -AND  PRINTED  PAGES  1  TO  2. 

-*8  PROC  PRI.NTTO  NEW  UNIT-20; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 

49  PROC  SORT  OUT=SCTR; 

50  BY  2  X  Y; 

NOTE:  DATA  SET  WORK.SCTR  HAS  3096  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  PROCEDURE  SORT  USED  0.76  SECONDS  .AND  6928K. 

51  PROC  SUMMARY; 

52  BY  2  X; 

53  VAR  Y; 

54  OUTPUT  OUT-OVLMN  MEAN-MEAN; 

NOTE:  THE  DATA  SET  WORK.OVLMN  HAS  48  OBSERVATIONS  .AND  5  VARIABLES. 

NOTE:  THE  PROCEDURE  SUMMARY  USED  0.56  SECONDS  AND  464K. 

55  data  sepin: 

56  SET  SCTR  OVL.MN ; 

5"  BY  2; 

NOTE:  DATA  SET  WORK. SEPIN  HAS  3144  OBSERVATIONS  .AND  10  VARIABLES . 

NOTE:  THE  DATA  STATEMENT  USED  0.67  SECONDS  AND  336K. 

58  PROC  PLOT  NOLEGEND  DATA-SEPIN; 

59  BY  2; 

60  PLOT  MEAN*X- ’ X '  Y*X-'.'  /  HAXIS— 3  TO  20  BY  1  VaXIS-50  TO  250  BY  25  OVERLAY 

61  ; 

62  TITLE  'EPINEPHRINE  SCATTER  DIAGRAM'; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.65  SECONDS  .AND  464K  AND  PRINTED  PAGES  3  TO  4 . 

63  PROC  PRINTTO  NEW  UNIT-21; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  .AND  336K. 

6**  PROC  PLOT  NOLEGEND  DATA-SEPIN; 

65  PLOT  MEaN-'X*  '  X '  /  HAXIS--3  TO  20  BY  1  VaXIS-50  TO  250  3Y  25; 

66  TITLE  'Mean  Epinephrine  Concentration  Versus  Time'; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.46  SECONDS  .AND  464K  .AND  PRINTED  PAGE  5. 

6'  PROC  PRINTTO  NEW  UNIT-22; 

63  TITLE  'CATECHOLAMINE  ANALYSIS:  Epinephrine’; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 

69  PROC  DATASETS; 

70 

LIST  OF  MEMBERS  BEFORE  UPDATE  OF  DIRECTORY. 

NAME  MEMTYPE  DBS  TRACKS  PPOT 

OVL.MN  DaTa  -3  1 


3  SASCR;  LOG'  CMS  SAS  5.16  VM/CMS  CMS  USER  QSECLSB 

TESTE  /DATA  3096  1 

TO  DELETE  SCTR : 

71  DELETE  OVLMN; 

LIST  OF  MEMBERS  AFTER  UPDATE  OF  DIRECTORY. 

NAME  MEMTYPE  OBS  TRACKS  PROT 

SEP  IN  / DATA  3144  1 

TESTE  DATA  3096  1 

NOTE:  THE  PROCEDURE  DATASETS  USED  0.12  SECONDS  AND  464K. 

"2  PROC  STEPWISE: 

'3  MODEL  Y  =  X  XSQR  Z  XZ  XSORZ  /  SLENTRY=0 . 095  SLST.AY-0.095  STEPWISE  MAXR ; 

NOTE:  THE  PROCEDURE  STEPWISE  USED  0.57  SECONDS  AND  464K  .AND  PRINTED  PAGES  6  TO  8 . . 


74  PROC  PRINTTO  NEW  UNIT-23; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 

75  PROC  REG; 

76  MODEL  Y  -  X  XSQR  /  PARTIAL; 

77  ID  CASE; 

NOTE:  ACOV  AND  SPEC  OPTION  ONLY  VALID  WITH  RAWDATA 

NOTE:  THE  PROCEDURE  REG  USED  1.76  SECONDS  AND  656K  .AND  PRINTED  PAGES  9  TO  12. 

78  PROC  PRINTTO  NEW  UNIT-24; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 

79  PROC  GLM; 

80  CLASS  X  Z; 

81  MODEL  Y  -  X  X*X  X*Z; 

NOTE:  THE  PROCEDURE  GLM  USED  3.10  SECONDS  AND  1040K  AND  PRINTED  PAGES  13  TO  14. 

82  PROC  PRINTTO  NEW  UNIT-25; 

83  * - 

84  * 

85  *  to  obtain  tables  listing  the  variance  inflation  factors, 

86  :f  influence  s  tat  i  s  t  i  cs  ,  and  tolerances,  the  following  SAC 

87  ”  statements  were  used  in  this  partition: 

88  * 

89  *  PROC  REG; 

90  MODEL  Y  =  X  XSQR  /  TOL  VIF  INFLUENCE; 

91  *  ID  CASE; 

92  "  OUTPUT  OUT-REPIN  P-PREDICT  R-RESID  STUDENT-STUDENT; 

93 

94  K - 


NOTE:  THE  PROCEDURE  PRINTTO  USED  0.04  SECONDS  AND  336K. 

--  ?RIC  REG; 

MODEL  Y  -  X  XSQR  /  I  SSI  SS2  STB  COVB  CORRB  SEQB  COLLIN 

COLLINOINT  ACOV  P  R  CLM; 

ID  CASE; 

lUTPUT  OUT-REPIN  P-PREDICT  R-RESID  STUDENT-STUDENT; 

' 'I  OaTa  SET  WORK. REPIN  HAS  3144  OBSERVATIONS  AND  13  VARIABLES. 

ICEDl'RE  REG  USED  6.93  SECONDS  AND  656K  AND  PRINTED  PAGES  15  TO  82. 
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4  SAS(R)  LOG  CMS  SAS  5.16  VM/CMS  CMS  USER  QSECLSB 

MOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 


101  PROC  PLOT  DATA” REP  IN ; 

102  PLOT  RES ID"X= ' " '  /  HAXIS— 3  TO  20  BY  1  V  AXIS— 150  TO  150  BY  25; 

103  PLOT  RES ID:''PREDICT= 1  * ’  /  HAXIS-115  TO  185  BY  5  V  AXIS— 150  TO  150  BY  25: 

104  PLOT  STUDENT*X« ’ * '  /  HAXIS— 3  TO  20  BY  1  V  AXIS  — 4  TO  5  BY  0.5; 

105  PLOT  STUDENT5tPREDICT= '  * '  /  HAXIS-115  TO  185  BY  5  V  AXIS  — 4'  TO  5  BY  0.5; 

106  TITLE  'EPINEPHRINE  RESIDUAL  PLOTS'; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.95  SECONDS  AND  464K  AND  PRINTED  PAGES  83  TO  86. 

107  PROC  PRINTTO  NEW  UNIT-27; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 

108  PROC  PLOT  DATA” RE? IN ; 

109  BY  2; 

110  PLOT  RES ID*CASE- ' * '  /  HAXIS-0  TO  65  BY  5  VAXIS— 150  TO  150  BY  25; 

111  PLOT  STUDENT*CASE= ' * '  /  HAXIS-0  TO  65  BY  5  V AXIS— 4  TO  5  BY  0.5; 

112  TITLE  ’EPINEPHRINE  RESIDUAL  PLOTS'; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.79  SECONDS  AND  464K  AND  PRINTED  PAGES  87  TO  90. 

113  PROC  PRINTTO  NEW  UNIT-28; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  336K. 


114  PROC  AUTOREG; 


115 

TITLE 

' Epinephrine 

Autoregressive  Models' 

1 

116 

MODEL 

Y 

«  X 

XSQR 

/ 

COEF  CQRRB  COVB  BACKSTEP; 

117 

MODEL 

Y 

-  X 

XSQR 

/ 

NLAG-1  COEF  CORRB 

COVB 

BACKSTEP ; 

118 

MODEL 

Y 

-  X 

XSQR 

/ 

NLaG-2  COEF  CORRB 

COVB 

BACKSTEP; 

119 

MODEL 

Y 

=  X 

XSQR 

/ 

NLAG-3  COEF  CORRB 

COVB 

BACKSTEP ; 

120 

MODEL 

Y 

-  X 

XSQR 

/ 

NLAG-4  COEF  CORRB 

COVB 

BACKSTEP ; 

NOTE: 

THE  PROCEDURE  AUTOREG 

USED  6.92  SECONDS 

AND 

464K  AND  PRINTED 

NOTE:  SAS  USED  6928K  MEMORY. 


NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 

PO  BOX  8000 

CARY,  N.C.  27511-8000 
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SAS(R)  LOG  CMS  SAS  5.16 


VM/ CMS  CMS  USER  QSECLSB 


NOTE:  COPYRIGHT  (C)  1984,1986  SaS  INSTITUTE  INC.,  CARY,  N.C.  27511,  U.S.A. 
NOTE:  CMS  SaS  RELEASE  5.16  AT  GEORGIA  INSTITUTE  OF  TECHNOLOGY  (03559001). 

NOTE:  CPU ID  VERSION  -  FF  SERIAL  -  012242  MODEL  -  4381  . 

NOTE:  SAS  OPTIONS  SPECIFIED  ARE: 

LEAVE-0 

1  DATA  TESTD; 

2  CMS  FILEDEF  X  DISK  DO PAM IN  DAT  Al ; 


3 

CMS 

FILEDEF 

20 

DISK 

DOPAMINO 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

<4 

CMS 

FILEDEF 

21 

DISK 

DOPAMIN1 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

5 

CMS' 

’FILEDEF 

22 

DISK 

DOPAMIN2 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

6 

CMS 

FILEDEF 

23 

DISK 

DOPAMIN3 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

7 

CMS 

FILEDEF 

24 

DISK 

D0PAMIN4 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

8 

CMS 

FILEDEF 

25 

DISK 

DO PAM I N5 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

9 

CMS 

FILEDEF 

26 

DISK 

DO PAM I N6 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

10 

CMS 

FILEDEF 

27 

DISK 

DOPAMIN7 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

11 

CMS 

FILEDEF 

28 

DISK 

DOPAMIN8 

LISTING 

Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

12 

ARRAY  WEEK 

{24} 

WKN3  WKN2  MISSN1  WKO- 

-WK20; 

13  KEEP  X  XSQR  Y  2  XZ  XSQRZ  CASE; 

14  INFILE  X; 

15  INPUT  CASE  1-3 


16 

WKN3 

5-7 

17 

WKN2 

9-11 

18 

WKO 

13-15 

19 

WK1 

17-19 

20 

WK2 

21-23 

21 

WK3 

25-27 

22 

WK4 

29-31 

23 

WK5 

33-35 

24 

WK6 

37-39 

25 

WK7 

41-43 

26 

WK8 

45-47 

27 

WK9 

49-51 

28 

WK10 

53-55 

29 

WK1 1 

57-59 

30 

WK12 

61-63 

31 

WK13 

65-67 

32 

WK1 4 

69-71 

33 

WK15 

73-75 

34 

WK16 

77-79 

35 

WK17 

81-83 

36 

WK18 

85-87 

37 

WK19 

89-91 

38 

WK20 

93-95 

39  ; 

<*0  MISSN1- . ; 

-1  IF  CASE  <  100  THEN  Z  -  0; 

42  IF  CASE  >-  100  THEN  Z  -  1; 

43  IF  Z-l  THEN  CASE-CASE- 100; 

44  do  I  -  1  TO  24; 

45  X  -  1-4;  XSQR  -  X*X;  XZ  -  X*Z;  XSQRZ  -  X*X*Z;  Y  -  WEEK  tl) ; OUTPUT; 

46  END; 


NOTE:  INFILE  X  IS  FILE  DC PAM IN  Da7  a! 
NOTE:  128  -INES  WERE  READ  FROM  INFILE  X. 
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SAS(R)  LOG  CHS  SAS  5.16 


VM/CMS  CMS  USER  QSECLSB 


NOTE:  DATA  SET  WORK. TESTD  HAS  3072  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.58  SECONDS  AND  252K. 

47  PROC  CONTENTS; 

NOTE:  THE  PROCEDURE  CONTENTS  USED  0.20  SECONDS  AND  316K  AND  PRINTED  PAGES  1  TO  2. 

48  PROC  PRINTTO  NEW  UNIT-20; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 

49  PROC  SORT  OUT-SCTR; 

50  SY  Z  X  Y; 

NOTE:  DATA  SET  WORK.SCTR  HAS  3072  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  PROCEDURE  SORT  USED  0.78  SECONDS  AND  6908K. 

51  PROC  SUMMARY; 

52  BY  Z  X; 

53  VAR  Y; 

54  OUTPUT  OUT-OVLMN  MEAN-MEAN; 

NOTE:  THE  DATA  SET  WCRK.OVLMN  HAS  48  OBSERVATIONS  AND  5  VARIABLES. 

NOTE:  THE  PROCEDURE  SUMMARY  USED  0.57  SECONDS  AND  444K. 

55  DATA  SDOPAMIN; 

56  SET  SCTR  OVL.MN; 

57  BY  Z; 

NOTE:  DATA  SET  WORK. SDOPAMIN  HAS  3120  OBSERVATIONS  AND  10  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.57  SECONDS  AND  316K. 

58  PROC  PLOT  NOLEGEND  DATA-SDOPAMIN; 

59  BY  Z; 

60  PLOT  MEAN*X- ' X '  Y*X-’.’  /  HAXIS— 3  TO  20  BY  1  VAXIS-0  TO  100  BY  10  OVERLAY; 

61  TITLE  'DOPAMINE  SCATTER  DIAGRAM’; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.66  SECONDS  AND  444K  AND  PRINTED  PAGES  3  TO  4. 

62  PROC  PRINTTO  NEW  UNIT-21; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 

63  PROC  PLOT  NOLEGEND  DaTA-SDOPAMIN; 

64  PLOT  MEAN*X-'X'  /  HAXIS— 3  TO  20  BY  1  VAXIS-0  TO  100 • BY  10; 

65  TITLE  'Mean  Dopamine  Concentration  Venus  Time'; 

NOTE:  THE  PROCEDURE  PLOT  USED  0.47  SECONDS  AND  444K  AND  PRINTED  PAGE  5. 

66  PROC  PRINTTO  NEW  UNIT-22; 

67  TITLE  'CATECHOLAMINE  ANALYSIS:  Dopamine'; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.03  SECONDS  AND  316K. 


68  PROC  DATASETS; 

69 


LIST  OF 

MEMBERS 

BEFORE  UPDATE  OF 

DIRECTORY. 

NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

OVLMN 

/DATA 

48 

1 

SCTR 

/DATA 

3072 

1 

SDOPAMIN 7 DaTa 

3120 

1 

TESTD 

/  Da  t  a 

3072 

1 
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69  DELETE  SCTR; 

70  DELETE  OVLMN; 

LIST  OF  MEMBERS  AFTER  UPDATE  OF  DIRECTORY. 


NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

SDOPAMIN/DATA 

3120 

1 

TESTD 

/DATA 

3072 

1 

NOTE:  THE  PROCEDURE  DATASETS  USED  0.12  SECONDS  AND  444K. 

71  PROC  STEPWISE; 

72  MODEL  Y  -  X  XSQR  Z  XZ  XSQRZ  /  SLENTRY-0.10  SLSTAY-0.10  STEPWISE  MAXR ; 
NOTE:  THE  PROCEDURE  STEPWISE  USED  0.60  SECONDS  AND  444K  AND  PRINTED  PAGES  6  TO  9 

73  PROC  PRINTTO  NEW  UNIT-23; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 

74  PROC  REG; 

75  MODEL  Y  -  X  XSQR  XZ  /  PARTIAL; 

76  '  ID  CASE; 

NOTE:  ACOV  AND  SPEC  OPTION  ONLY  VALID  WITH  RAWDATA 

NOTE:  THE  PROCEDURE  REG  USED  1.64  SECONDS  AND  636K  AND  PRINTED  PAGES  10  TO  14. 

77  PROC  PRINTTO  NEW  UNIT-24; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 

78  PROC  GLM; 

79  CLASS  X  Z; 

80  MODEL  Y  -  X  X*X  X*Z; 

NOTE:  THE  PROCEDURE  GLM  USED  3.20  SECONDS  AND  1020K  AND  PRINTED  PAGES  15  TO  16. 

81  PROC  PRINTTO  NEW  UNIT-25; 

82  * - 

83  * 

84  *  to  obtain  tables  listing  the  variance  inflation  factors, 

85  *  influence  statistics,  and  tolerances,  the  following  SAS 

86  *  statements  were  used  in  this  partition: 

87  * 

88  *  PROC  REG; 

89  *  MODEL  Y  -  X  XSQR  XZ  /  TOL  VIF  INFLUENCE; 

90  *  ID  CASE; 

91  *  OUTPUT  OUT-RDOPAMIN  P-PREDICT  R-RESID  STUDENT-STUDENT; 

92  * 

93  - 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.04  SECONDS  AND  316K. 

94  PROC  REG; 

95  MODEL  Y  -  X  XSQR  XZ  /  I  SSI  SS2  STB  COVB  CORRB  SEQB  COLLIN 

96  COLLINOINT  ACOV  P  R  CLM; 

97  ID  CASE; 

98  OUTPUT  OUT-RDOPAMIN  P-PREDICT  R-RESID  STUDENT-STUDENT; 

NOTE:  THE  DATA  SET  WORK.RDOPAMIN  HAS  3120  OBSERVATIONS  AND  13  VARIABLES. 

NOTE:  THE  PROCEDURE  REG  USED  7.33  SECONDS  AND  636K  AND  PRINTED  PAGES  17  TO  83. 

99  PROC  PRINTTO  NEW  UNIT-26: 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 
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100  PROC  PLOT  DATA-RDOPAMIN; 

101  PLOT  RESID,fX*  ’  * '  /  HAXIS— 3  TO  20  BY  1  VAXIS— 125  TO  125  BY  25; 

102  PLOT  RESID*PREDICT-’*'  /  HAXIS-15  TO  65  BY  5  VAXIS— 125  TO  125  BY  25; 

103  PLOT  STUDENT*X- '  * '  /  HAXIS— 3  TO  20  BY  1  VAXIS— 2  TO  6  BY  0.5; 

104  PLOT  STUDENT*PREDICT- ' * ’  /  HAXIS-15  TO  65  BY  5  VAXIS— 2  TO  6  BY  0.5; 

105  TITLE  ' DOPAMINE  RESIDUAL  PLOTS’; 

MOTE:  THE  PROCEDURE  PLOT  USED  0.96  SECONDS  AND  444K  AND  PRINTED  PAGES  84  TO  87. 

106  PROC  PRINTTO  NEW  UNIT-27; 

MOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 

107  PROC  PLOT  DATA-RDOPAMIN; 

108  BY  2; 

109  PLOT  RES ID*CASE“ ’ * '  /  HAXIS-0  TO  65  BY  5  VAXIS— 125  TO  125  BY  25; 

110  PLOT  STUDENT "CASE- ’ * 1  /  HaXIS-O  TO  65  BY  5  VAXIS— 2  TO  6  BY  0.5; 

111  TITLE  'DOPAMINE  RESIDUAL  PLOTS'; 

NOTE:  THE  PROCEDURE  PLOT  USED'0.79  SECONDS  AND  444K  AND  PRINTED  PAGES  88  TO  91. 

112  PROC  PRINTTO  NEW  UNIT-28; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  316K. 

113  PROC  AUTOREG; 

114  TITLE  'Dopamine  Autoregressive  Models’; 

115  MODEL  Y  -  X  XSQR  XZ  /  CO EF  CORRB  COV B  BACKSTEP; 

116  MODEL  Y  -  X  XSQR  XZ  /  NLAG-1  COEF  CORRB  COVB  BACKSTEP; 

117  MODEL  Y  -  X  XSQR  XZ  /  NLAG-2  COEF  CORRB  COVB  BACKSTEP; 

118  MODEL  Y  -  X  XSQR  XZ  /  NLAG-3  COEF  CORRB  COVB  BACKSTEP; 

119  MODEL  Y  -  X  XSQR  XZ  /  NLAG-4  COEF  CORRB  COVB  BACKSTEP; 

NOTE:  THE  PROCEDURE  AUTOREG  USED  6.82  SECONDS  AND  444K  AND  PRINTED  PAGES  92  TO  104. 
NOTE:  SAS  USED  6908K  MEMORY. 

NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 

PO  BOX  8000 

CARY,  N.C.  27511-8000 
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SUMMARY  OF  S1EPW1SE  REGRESSION  PROCEDURE  FUR  OEPENUENI  VARIAUIE 
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GENERAL  LINEAR  MODELS  PROCEOURE 
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^this  term  contains  both  sum-of-squares  pure  error  and 
1.38  sum-of-squares  lack-of-fit. 
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